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Abstract 

Direct numerical simulations of the interaction of a premixed flame with driven, subsonic, homogeneous, isotropic, 
Kolmogorov-type turbulence in an unconfined system are used to study the mechanisms determining the turbulent 
flame speed, St, in the thin reaction zone regime. High intensity turbulence is considered with the r.m.s. velocity 35 
times the laminar flame speed, S l, resulting in the Damkohler number Da - 0.05. The simulations were performed 
with Athena-RFX, a massively parallel, fully compressible, high-order, dimensionally unsplit, reactive-flow code. 
A simplified reaction-diffusion model, based on the one-step Arrhenius kinetics, represents a stoichiometric H2-air 
mixture under the assumption of the Lewis number Le - \. Global properties and the internal structure of the 
flame were analyzed in an earlier paper, which showed that this system represents turbulent combustion in the thin 
reaction zone regime. This paper demonstrates that: (1) The flame brush has a complex internal structure, in which 
the isosurfaces of higher fuel mass fractions are folded on progressively smaller scales. (2) Global properties of 
the turbulent flame are best represented by the structure of the region of peak reaction rate, which defines the flame 
surface. (3) In the thin reaction zone regime, is predominantly determined by the increase of the flame surface 
area, Aj, caused by turbulence. (4) The observed increase of St relative to S l exceeds the corresponding increase of 
At relative to the surface area of the planar laminar flame, on average, by ^ 14%, varying from only a few percent to as 
high as ^30%. (5) This exaggerated response is the result of tight flame packing by turbulence, which causes frequent 
flame collisions and formation of regions of high flame curvature > 1 jSi, or "cusps," where Sl is the thermal width 
of the laminar flame. (6) The local flame speed in the cusps substantially exceeds its laminar value, which results in a 
disproportionately large contribution of cusps to 5 7 compared with the flame surface area in them. (7) A criterion is 
established for transition to the regime significantly influenced by cusp formation. In particular, at Karlovitz numbers 
Ka > 20, flame collisions provide an important mechanism controlling St, in addition to the increase of At by 
large-scale motions and the potential enhancement of diffusive transport by small-scale turbulence. 
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1. Introduction 

One of the fundamental questions of turbulent combustion research concerns our ability to understand and predict 
the rate of energy release, or equivalently the burning speed, of a turbulent flame. This can be achieved if two 
key aspects of the combustion process can be determined: (1) its local properties, namely the local speed of flame 
propagation, and (2) its global characteristics, i.e., the overall structure of the turbulent flame that connects the local 
burning velocity with the turbulent flame speed. In general, such local and global characteristics are not universal for 
the combustion process. They can vary substantially both in space and time due to the unsteadiness, inhomogeneity, 
and anisotropy of the turbulent flow associated with the system geometry, presence of walls and boundaries, change 
in the flow conditions, etc. 

Our present understanding of turbulent combustion is largely based on the concept originally proposed almost 
70 years ago by Damkohler |1] (also see reviews by |<2|,|3i|4|,|51]). This concept is formulated for a turbulent flame 
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that is, on average, a planar, quasi-one-dimensional structure propagating in a steady, homogeneous, and isotropic 
background turbulent flow. According to Damkohler's suggestion, two qualitatively distinct regimes of turbulence- 
flame interaction can be identified, and they are typically referred to as "large-scale" and "small-scale" turbulence 

m. 

In the first regime, the overall structure of the turbulent flame is determined by turbulent motions on large scales 
greater than the flame width, 6l- These motions stretch and fold the flame, thereby increasing its surface area, Aj- At 
the same time, they do not affect the local burning velocity, S„, with which the flame propagates normal to its surface 
at each point and which remains equal to the speed, S l, of the unperturbed laminar flame. Consequently, the speed, 
St, of the turbulent flame grows proportionally to the increase of the flame surface area, i.e., 5 7- = 5 l{At /Al), where 
Al is the surface area of the planar laminar flame. The limitations of this assumption were known practically at the 
time when this model was proposed |6]. In particular, when the Lewis number Le + 1, 5„ can vary with the curvature 
of the flame. As a result, the stretching and folding action of scales X> dL ends up affecting Sn- The relation between 
Aj and Sj, however, can be generalized by incorporating the eff'ects of flame strain and curvature in terms of the 
stretch factor / ([jst], also see Q for the review of the theory of flame stretch). 

In the second regime, large-scale wrinkling of the flame by turbulence is absent. Instead, small-scale motions, 
which are energetic enough, penetrate and disrupt the internal flame structure. Thus molecular diffusivity and thermal 
conduction become enhanced, or completely dominated, by the turbulent transport associated with such small-scale 
turbulence. As a result, in this regime, the turbulent flame speed is suggested to be HI] 
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where D and Dj are, respectively, the molecular and the effective turbulent diffusivities. 

These two regimes, in eff'ect, represent two limits of low and high turbulent intensity. In the first case, large-scale 
turbulence must be slow enough so that the Kolmogorov scale rj > 6l, thus causing the flow to remain laminar on 
scales A < Sl. This combustion regime corresponds to the "wrinkled" and "corrugated flamelets." In the opposite 
limit, the small-scale turbulence regime is realized when turbulent motions on all scales in the system are sufficiently 
fast to form an effectively one-dimensional (ID) structure analogous to the planar laminar flame, but which has a much 
larger width and speed. This second regime is often associated with either a "thin reaction zone" flame I5i] or a "broken 
reaction zone" ("distributed") flame UHHIifl]- As a consequence of considering these two limits, description of the 
turbulent flame becomes significantly simplified. The large-scale regime assumes very simple local flame properties, 
i.e., those of the laminar flame, but a complex global structure of a highly convolved flame surface. The small-scale 
regime, on the other hand, trades the simplicity of the global structure of a planar flame for the complexity of its 
internal structure, which is pronouncedly distinct from that of the laminar flame. Therefore, each of these regimes 
has only one dynamically important quantity that completely determines St, namely At for the large-scale turbulence 
and Dt for the small-scale turbulence. 

This picture then naturally raises the following question. How does the transition occur between these two limiting 
regimes of combustion? At some intermediate turbulent intensities, the flow will cease to be laminar on scales A < di. 
As a result, turbulence will begin to penetrate the flame, disrupting it and causing its internal structure to differ from 
that of the laminar flame. At the same time, the flame can remain sufficiently thin to form, under the action of high- 
intensity large-scale motions, a highly convolved and tightly packed configuration. The picture outlined above would 
then imply a gradual shift with increasing turbulent intensity from the large-scale turbulence being the dominant 
process controlling 5 7 to the small-scale one. Over a substantial range of turbulent intensities, however, it must be 
their combined action that determines the magnitude of St- 

Of particular interest, in this respect, is the "thin reaction zone" regime, in which only the preheat zone is disrupted 
by turbulence. The presence of turbulent motions on scales A < 6l suggests that the molecular and thermal transport 
between the reaction and preheat zones are affected by turbulence. This should cause S„ to differ from St, which 
is characteristic of the small-scale turbulence regime. On the other hand, if the structure of the reaction zone is the 
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same as in the laminar case, as occurs at the large-scale boundary of this regime, then S„ must equal S l. Indeed, 
Sn = (1/po) / pYd^, where po is the fuel density, Y is the reaction rate, and ( is the coordinate normal to the flame. 
Since Y Q only in the reaction zone where it does not differ from the laminar flame, it follows that S„ = S l. What 
does then determine St in the thin reaction zone regime? Is it folding and stretching of the flame by turbulence, 
implying the validity of eq. ([TJ, or the microscopic turbulent transport represented by eq. or both? How does the 
balance between the two processes change with turbulent intensity? Finally, is the eff'ect of high-speed turbulence on 
the flame limited to the two processes discussed above or do other mechanisms also play a role? 

Before these questions can be addressed, several conceptual difficulties must first be resolved. Most importantly, 
the definition of the flame surface and its area must be revisited. This question, which has a seemingly simple answer 
at low turbulent intensities, becomes quite nontrivial when the flow is no longer laminar on scales A < 6l- In this 
case, different regions of the flame are affected differently by turbulence and, thus, the isosurfaces of temperature or 
fuel mass fraction are not necessarily parallel. Furthermore, without an accurate definition of the flame surface, it is 
not clear what the local burning speed of the flame is, since this quantity typically describes flame propagation normal 
to its surface and determining its local value requires integration of the distribution of Y normal to the flame surface, 
as was discussed above. 

In recent years, the advent of direct numerical simulations (DNS) has provided a tool permitting the mechanisms 
that determine 5 r to be directly investigated. DNS have been used by a number of groups to study the dynamics and 
properties of turbulent flames |9]-| 11] (also see a review by li2Zl ). With the exception of the work by Aspden et al. 
i9l ll0ll that analyzed the fully distributed thermonuclear flame, these studies generally assume relatively low turbulent 
intensities. Several authors also did explicitly consider the relation between St and At- Bell et al. ll20ll analyzed 
the three-dimensional (3D) DNS of a statistically planar turbulent methane flame and determined that / deviates 
from unity by < 10%. A similar estimate was obtained by Khokhlov ll23ll in 3D simulations of Rayleigh-Taylor 
driven thermonuclear flames in degenerate matter Both of these studies, however, considered reactive mixtures in 
which Le + I and, therefore, / would not be expected to be exactly equal to one. In this context, the result of 
1 2311 is particularly interesting since degenerate matter is characterized by Le » 1 and typically in excess of 10^. 
Furthermore, both Bell et al. lEoll and Khokhlov considered turbulent intensities which were too low to allow 
small-scale motions to penetrate efficiently the interior of the flame. For instance, in |20], the characteristic turbulent 
velocity was ^4.35'z,, and in lEsll . it reached values < 125'l. Consequently, no effects of the microscopic turbulent 
transport would be expected in these cases. 

In this respect, the work of Hawkes & Chen 112 ill is of particular interest as it considered the validity of eq. ([1]) in the 
thin reaction zone regime. They analyzed simulations of lean, statistically flat turbulent methane flames at turbulent 
intensities < 28.55 z, and found values of / within a few percent of unity, with the exception of the methane-hydrogen 
flame for which / = 1.18 was determined. Interpretation of the implications of these results for the thin reaction zone 
regime, however, is complicated by the two-dimensional (2D) nature of their simulations. Moreover, the turbulent 
integral scale was < Si, which substantially suppressed flame wrinkling and prevented the development of a highly 
convolved flame which would typically be expected to form in a realistic setting. 

The studies performed to-date, thus, generafly show that, at low turbulent intensities, the increase in 5 j- is indeed 
almost completely determined by the growth of At- At the same time, they leave open the questions discussed above 
regarding the mechanisms controlling in the thin reaction zone regime, let alone at higher turbulent intensities. 

The objective of this paper is to begin addressing these questions systematically by first considering the interaction 
of a statistically planar, turbulent flame with the steady, homogeneous, isotropic turbulence. We assume sufficiently 
high turbulent intensity to represent the regime that is borderline between thin and broken reaction zones, according 
to the traditional combustion regime diagrams 24, 25] (cf. Fig. 2 in [11]). In other words, we consider the highest- 
intensity turbulence, which has been hypothesized to allow for the existence of a flame with the internal structure of 
the reaction zone that is essentially not affected by turbulent transport. The turbulent rm.s. velocity of the flow in 
cold fuel is Urms ~ 355'l, leading to the Damkohler number Da = 0.05 and Gibson scale Lc ~ 3 x 10 '*5z,. Here we 
focus on the Le - I situation to exclude thermodiffiisive effects as a potential source of variations of the local flame 
speed. Developing rigorous understanding of this simpler case will then serve as a starting point for the study of more 
realistic situations represented by more complex chemistry. 

This paper continues the analysis of the simulations first presented in fv\^- The primary focus of flT] was a de- 
tailed study of flame properties and evolution in the presence of such high-speed turbulence. In particular, a method 



3 



was presented which allowed the direct determination of the internal structure of the flame based on its actual 3D 
configuration inside the flame brush. The analysis showed that the preheat zone is broadened by turbulence while the 
reaction-zone structure remains virtually identical to that of the planar laminar flame. Therefore, this study demon- 
strated that, even under the action of such intense turbulence, the system can be classified as being in the thin reaction 
zone regime. 

Our starting point in this work is the question of the definition of the flame surface area in the high-speed turbulent 
flow. Answering this will then allow us to determine whether a relation between St and At, similar to eq. ([T]i, can 
indeed be established and, thus, which processes control 5 j- in the thin reaction zone regime. 



2. Numerical method and simulations performed 



2.1. Physical and numerical models 

Here we summarize the physical model, the numerical method used, and key aspects of the simulation setup. A 
more detailed discussion can be found in ifHIl . 

We solve the system of unsteady, compressible, reactive flow equations. 



(3) 
(4) 
(5) 
(6) 



Here p is the mass density, u is the velocity, E is the energy density, P is the pressure, Y is the mass fraction of the 
reactants, q is the chemical energy release, and Y is the reaction source term. The coefficients of thermal conduction, 
TC, and molecular diffusion, D, are 

7^" '7C 7^" 
D = Do — , — = ^0 — , (7) 
P pCp p 

where Do, a-o, and n are constants, and Cp - jR/Miy - 1) is the specific heat at constant pressure. The equation of 
state is that of an ideal gas. Chemical reactions are modeled using the first-order Arrhenius kinetics 
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where B is the pre-exponential factor, Q is the activation energy, and R is the universal gas constant. 

Table |2] summarizes the parameters of the physical model used as well as the resulting properties of the planar 
laminar flame Ol ill . These parameters are based on a simplified reaction-diffusion model designed to represent the 
stoichiometric H2-air mixture |26]. We refer to 111, 2^ for the detailed discussion of the limitations of applicability 
of this model for the description of the actual Ha-air mixture. 

Eqs. (O-® are solved using the code Athena-RFX flT] - the reactive flow extension of the magnetohydrody- 
namic code Athena (rfJlS]. Athena-RFX is a fixed-grid, massively parallel, finite-volume, fully conservative code. It 
implements a variant [28.1 of the fully unsplit corner-transport upwind (CTU) algorithm |29] and its 3D extension pre- 
sented in isoll . in conjunction with the PPM spatial reconstruction llsill and the approximate nonlinear HLLC Riemann 
solver Overall, the code achieves 3''''-order accuracy in space and 2'"'-order accuracy in time. A detailed description 
and an extensive suite of tests of the hydrodynamic integration algorithm can be found in 11271 12811 . Implementation 
of the reactive-difiiisive extensions in Athena-RFX, along with the results of tests including convergence studies, is 
discussed in detail in ll lll32ll . 

The simulations presented here model flame interaction with steady, homogeneous, isotropic turbulence, described 
by the classical Koknogorov theory ll33ll . Turbulence driving is implemented using a spectral method | U , 32 ] similar 



to the one used in |34, 35[]. In this method, velocity perturbations 5u(k) are initialized in Fourier space, and each 
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component 5Mj(k) is an independent realization of a Gaussian random field superimposed with the desired energy 
injection spectrum. Subsequently, nonsolenoidal components of 5u(k) are removed to ensure that V ■ 5u(x) - 0. An 
inverse Fourier transform of 5u(k) gives 5u(x), the velocity perturbation field in the physical space. The Su{x) is 
normalized to provide a constant rate e of kinetic-energy injection, and the total momentum in the perturbation field is 
subtracted from <5u(x) to ensure that no net momentum is added to the domain, i.e., j pSu - 0. The resulting velocity 
perturbations are added to the flow field u(x) at every time step, and the overall perturbation pattern is regenerated at 
every time interval Atvp ~ lOAx/c,, where Cj is the maximum sound speed in the domain and Ax is the computational 
cell size. 

Energy is injected only at the scale of the domain width, L, to obtain the Kolmogorov-type spectrum. The resulting 
turbulence is statistically steady, isotropic, and homogeneous with the inertial range of the energy cascade extending 
all the way to the energy injection scale (Fig.[TJ. Moreover, since the velocity perturbation field is divergence-free, no 
compressions or rarefactions are artificially induced as a result of driving. 

Eqs. (Is)-® do not explicitly include physical viscosity, but instead rely on numerical viscosity to provide the 
kinetic-energy dissipation. By systematically varying the numerical resolution in this approach, the energy spectrum 
can be extended to progressively smaller scales A «; (5l,o while maintaining the spectrum constant on larger scales. 
This is illustrated in Fig.[Tl which shows instantaneous kinetic energy spectra in the cold nonreactive flow immediately 
before the flame is initiated in the computational domain in three simulations with progressively increasing resolution. 
Such approach allows us to vary only the intensity of the small-scale turbulent motions, and thus to investigate their 
effects by dififerentiating them from the efi'ects of large scales. Furthermore, it shows the range of scales that must 
be resolved in a numerical simulation in order to capture accurately the evolution of the turbulent flame. This is 
particularly important in the context of high-speed turbulent reactive flows, in which it is typically impossible to 
resolve the Kolmogorov scale. The analysis presented in 111 ill demonstrated that scales A ^ ^^.o have virtually no 
effect on the properties of the turbulent flame with the exception of the degree of flame surface wrinkling on the fuel 
side. This showed that the evolution of the turbulent flame can be accurately reproduced without the need to resolve 
scales A «; ^^.o. We refer to 111 ill for a detailed discussion of this issue and, in particular, for a discussion of the 
applicability of the results to the actual stoichiometric H2-air mixture. 



2.2. Summary of simulations 

Key parameters of the simulations discussed here are summarized in Table [3] The main difference between the 
three calculations is the resolution, which progressively increases from Ax = 5l,o/8 in SI to Ax = <5l,o/32 in S3. The 
domain in all cases was initialized with uniform density po and temperature To (see Table|2]i. In S 1 and S2, initial fluid 
velocities were set to 0. In contrast, the velocity field in S3 was initialized with the ideal energy spectrum &(k) oc k^^/^ 
extending from the energy injection scale L to the numerical Kolmogorov scale 77 = 2Ax. This initial spectrum was 
normalized to ensure that at f = the total kinetic energy in the domain was equal to its predicted steady-state value 
ll32ll . In SI and S2, the flow field was allowed to evolve for the time tig„ - ItT^d, and in S3 for the time f,^„ - Ir^d, to 
develop the steady-state turbulent flow field. The corresponding instantaneous spectral density of the specific kinetic 
energy immediately prior to f,^„ in SI -S3 is shown in Fig. [1] At the time f,g„, a planar laminar flame with its front 
parallel to the x-y plan^H was initialized in the domain with the values of p, T, and Y based on the exact laminar flame 
solution. The velocity field was not modified, thereby preserving its structure that developed during the equilibration 
stage. 

Prior to ?,;,„, all domain boundaries were periodic. At tig,,, boundary conditions (BCs) along the left and right 
z-boundaries were switched to zero-order extrapolation in order to prevent pressure build-up inside the domain and 
to accommodate the resulting global fluid flow. Reflections of the acoustic perturbations caused by such BCs can be 
neglected in comparison with the acoustic noise generated by the turbulent flow field itself, and overall such choice of 
the BCs was found not to introduce any unphysical artifacts and thus not to aff'ect the solution accuracy. Furthermore, 
any potential effect of the outflow BCs was minimized by the large length-to- width ratio of the computational domain, 
which allowed the flame brush to remain sufficiently far from the boundaries at all times II ill . Hereafter, for simplicity 
we refer to the moment of ignition as f = 0. 



The longest dimension of the domain is along the z-axis. 
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Energy is injected into the domain with the constant rate e per unit volume for the total duration of the simulations 
to provide steady driving of the turbulent flow. The value of e was chosen to produce a high-intensity turbulence field 
that, however, was weak enough to minimize the probability of creating weak transonic shocklets arising from the 
intermittency in turbulent flow. The characteristic turbulent velocities and the turbulent integral scale of the steady- 
state flow before the moment of ignition are listed in Table [3] Corresponding values of the Damkohler number and 
the Gibson scale, given in the table, show that Da «; 1 and Lc 6l,o- Prior analysis of the resulting structure of 
the turbulent flame demonstrated that, despite the high turbulent intensity, small-scale turbulence fails to penetrate the 
internal structure of the reaction zone, and the system evolves in the thin reaction zone regime 

The following two key global characteristics of the turbulent flame are used extensively in this work. The width 
of the turbulent flame brush is defined as 

~ Zl,max ~ Z0,mm7 (9) 



where zo^min and zi^max are defined as 



- max(z) : Y(x,y,z) < 0.05 V (x,y,z < Zojnin), 
Zi,max = min(z) : Y(x,y,z) > 0.95 V (x,y,z > Zi,max)- 



(10) 



In other words, zo,min marks the rightmost x-y-plane, to the left of which is pure product, while zi^ma.x marks the 
leftmost x-y-plane, to the right of which is pure fuel. This is illustrated in Fig.|2] The turbulent flame speed is defined 
as 

where niR is the total rate of fuel consumption inside the flame brush, i.e., the total mass of reactants which is trans- 
formed to product per unit time. The detailed discussion of this choice of the definition of St can be found in 1 11 ] . 

Kinetic energy dissipation in the turbulent flow in the domain causes gradual heating of the fuel. The resulting 
relative increase in the internal energy, E,, and temperature of the nonreactive flow in one eddy turnover time, Ted, can 



be accurately estimated as 01 iL 13211 



Er — Et n T — Ta o 

= = DKrir - i)m4. (12) 

Here the constant Dk = 0.5 and £,,0 is the thermal energy in the domain at f = 0. Based on values of -y from Table|2]and 
Mof from Tabled T increases by ^ 1.8 K over Ted- Consequently, by the end of the simulation, the fuel temperature 
in the domain rises from To = 293 K to 320-330 k| As a result of such change in the thermodynamic state of 
fuel, the corresponding laminar flame properties also vary with time. In particular, this concerns the instantaneous 
thermal width, 61(1), and speed, S lit), of the laminar flame, which will deviate from their initial values 61^ and S'l.o 
corresponding to the fuel temperature To and pressure Pq (Table|2]i. 

In order to determine 6l and 5 l in simulations S 1 -S3, we found at each time Taverit) and Paverit) averaged over aU 
cells located between the planes z - Zi.max + L and z - Zi.max + 2L. This provided the average current thermodynamic 
state of fuel entering the flame brush. Based on Taver and Paver, an exact laminar flame solution was constructed, which 
gave the instantaneous values of Siit) and S lit). This analysis demonstrated that in the simulations, 6l decreases by 
< 15%. The S L changes more substantially increasing by < 25%. Thus, in order to consider accurately the nature of 
the turbulent flame speed in the high-speed turbulent reactive flows, it is important to account for the increase in the 
local burning speed of the flame due to the turbulent fuel heating. Consequently, in our further analysis, we relate S t 
to S L rather than S 1,0- 



Detailed discussion of the turbulent flame evolution in the simulations is given in 111 111 . Here, for completeness, 
we reproduce in Fig. |3]for all three calculations the time histories of 6t, normaUzed by 6l, as well as S t, normalized 
by5L(cf. Fig.4in[ll]). 

Finally, Table|4]lists the time-averaged values of both St/Sl andSr/S l [11] as well as of other key characteristics 
of the turbulent flame discussed below. Table |4] also shows the corresponding order of self-convergence for each 
quantity listed. Since the computational cell size decreases progressively by a factor of 2 for each simulation, the 



^Temperature increase occurs, on average, at constant density. 
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order of self-convergence (9(0) of a variable (p is defined as 



0(0) = log2 



Ifel -0531 
J0S2 -0S3l/ 



(13) 



3. Surface area and density of the fuel mass fraction isosurfaces 

3.1. Isosurface area 

Fig- llt>b shows the evolution of the surface areas Aq.oi and A0.99 of the isosurfaces of the fuel mass fraction 
Y - 0.01 and Y - 0.99, respectively. All isosurfaces are constructed using the "marching cubes" algorithm. These 
isosurfaces represent the outer boundaries that separate the flame brush from pure product and pure fuel. All values are 
normalized by the domain cross-section corresponding to the surface area of the planar laminar flame unperturbed 
by turbulence. In principle, the normalization should be performed over some average surface area of the flame brush, 
thus reflecting its overall large-scale shape. With suflicient accuracy, however, the turbulent flame in the system 
considered here can be viewed, on average, as a planar propagating front, and thus we normalize over its area L^. 

Fig. 3^ shows that Ao.oi exhibits substantial variability, similar to 6t and St (cf. Fig.O. There is a clear corre- 
spondence between the peaks and troughs of Ao.oi and St and, to a lesser extent, ^7-. Moreover, in a manner similar 
to 6t and St, Aqoi varies less with increasing resolution, and it appears on average to have converged. 

Such correlation between Ao.oi and St reflects the evolutionary cycle of the flame brush, discussed in fT\\\. In the 
turbulent flame, which on average is in a steady state, the influx of fresh fuel and the rate of its consumption never 
perfectly balance each other Instead, periodically either one or the other process dominates. In particular, an increase 
in the flame surface area inside the brush, and the associated increase in 6t, leads to a higher St- This causes rapid 
consumption of fuel inside the brush, decreasing both the flame surface area and the overall width of the flame brush. 
The result is a slower, less convolved flame, which is thus more susceptible to the action of turbulence. This leads to 
increased folding and stretching of the flame, and the cycle repeats. 

The behavior of A0.99 is pronouncedly diff'erent. There appears to be no correlation between variations in A0.99 
and either St ot St- Moreover, with increasing resolution, the magnitude of the variations increases, and there is no 
evidence of convergence of the growing average values of A0.99. We will consider the correlation between Ay and 
both St and 6t in a more quantitative form in §|4l 

These conclusions regarding the change of Aq.oi and A0.99 with resolution are supported by Fig.|4^, which shows 
the full normalized time-averaged distributions Ay/L^ for all values of Y. In simulation SI, Ao.oi and A0.99 are nearly 
equal. With increasing resolution, however, they diverge to the point that, in S3, there is almost a factor of 2 difference 
between them. Values of Ao.oi steadily decrease with resolution and they indeed demonstrate convergence. At the 
same time, A0.99 increases with resolution, and the diff'erence between S2 and S3 is only marginally smaller than 
between SI and S2. 

This behavior is part of a broader qualitative change in the overall shape of the distribution of Ay that occurs 
around Y ^ 0.5, i.e., at the boundary of the reaction and preheat zones. The lower resolution of SI suppresses small- 
scale turbulent motions, which, in turn, causes less wrinkling of the isosurfaces in the preheat zone. As a result, the 
flame-brush surface appears similar both on the product and fuel sides (cf. Fig. 3 in 11 111 ). With higher resolution, 
however, isosurface wrinkling becomes more pronounced with increasing distance from the reaction zone, and the 
distribution of Ay develops a distinctive inverted-S shape, as seen in Fig. |4^. The consequence of this was observed 
in 111 111 , which reported a much more highly convolved flame-brush surface on the fuel side in calculations S2 and S3. 
We will discuss the role of small-scale turbulence in further detail in § 13.41 

Profiles of Ay in S2 and S3 are closer to each other in the reaction zone than in the preheat zone (Fig.|4^). Table|4] 
shows that in the region of peak reaction rate, i.e., at F « 0.15, Ay/L^ exhibits the 3'^''-order convergence. This 
is substantially faster than what would be expected for an overall 2'"'-order-accurate code. Moreover, in all three 
calculations. Ay varies the slowest inside the reaction zone, i.e., for Y ^ 0.15 - 0.6. This shows that isosurfaces 
of Y, on average, follow each other relatively closely within the reaction zone, which thus can be viewed as being 
uniformly folded and stretched by turbulence as a coherent structure. This is consistent with the very low variability 
in the reaction zone of the instantaneous profiles of Y and T, which represent the internal flame structure fill] . 
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3.2. Isosurface density 

The analysis of the isosurface area distribution does not show how the flame is organized inside the flame brush 
and, in particular, how tightly it is packed. More specifically, the question arises whether variations in the isosurface 
area seen in Fig. |4^,b can be accounted for only by the change in the turbulent flame thickness, or whether they are 
also related to how the flame is folded inside the flame brush. 

Consider the surface density of the fuel mass fraction isosurfaces, defined as the isosurface area for a given Y 
normalized by the efi'ective flame-brush volume, 

(14) 



A detailed discussion of this definition of "Zy and, in particular, of the choice of the uniform normalization by the full 
flame-brush volume is given in |Appendix A| 

The evolution of l.om and S0.99 is shown in Fig.|4j;,d. The overall trends are similar to those exhibited by Aq.oi and 
A0.99. In particular, 1,y is lower on the product side of the flame. Similar to Ay, Sy rapidly converges with resolution 
at lower values of Y and shows virtually no convergence in the preheat zone. Most important, there is substantial 
variability in Sy, especially on the fuel side, where Sy changes by more than a factor of 3 in the highest resolution 
case. This tighter packing of the isosurfaces, which occurs periodically in the course of system evolution, suggests 
that the increase in Ay is associated not only with the increase of 6t but also of Zy. 

Fig. Hf shows the time-averaged distributions of Ey in simulations S 1-S3. Comparison with Fig. 15 demonstrates 
that both AyIL? and Sy show remarkably similar behavior The time-averaged values of Sy are also not constant, 
but, instead, for sufficiently high numerical resolution, they progressively increase through the flame interior from its 
product side to the fuel side. The resulting inverted-S shape of the profiles is very similar to that of Ay with both Ay 
and 2y at the two extreme values of F (F = 0.01 and 0.99) differing in S3 by almost a factor of two. Furthermore, 
similar to Ay, Sy varies the least in the interior of the reaction zone, i.e., for F » 0.15 - 0.6, while outside this region it 
shows a strong dependence on Y. Just outside the reaction zone, individual profiles begin to diverge from each other, 
with the variation between them becoming progressively larger for higher values of Y. 

In the reaction zone, values of 2y for all three calculations are much closer than those of Ay. In fact. Table |4] 
shows that Ey exhibits 4'''-order convergence compared with 3'''-order convergence for Ay. As a result, values of Ey 
are virtually identical for Y ^ 0.01 - 0.4 in S2 and S3. On the other hand, in the preheat zone, the profiles of Ey 
diverge more than in the case of Ay. This shows that the effects of small-scale turbulent motions on the reaction and 
the preheat zones are, in fact, even more disparate in terms of their ability to provide tight folding of the isosurfaces, 
compared with their ability to increase the isosurface area. 

Flame surface density is a quantity used both in experimental and theoretical combustion research. In particular, 
E,„fli is defined as the surface density in the region of the flame brush with mean reactedness c = 0.5, where c - 
(T — Tq)I(Tp - To) and Tp is the post-flame temperature. It can then be shown that 

X,„.. . (15) 

where A^^ is the flame surface area based on c = 0.5, and A^ represents the average area of the flame brush. In a 
system described by the one-step Arrhenius kinetics, c - 0.5 corresponds to F = 0.5. Therefore, A^^ = A0.5. At 
the same time, in a planar brush A^ = L^, as was discussed in § 13.11 Thus, I-max is equivalent to E0.5 as given by 
eq. (fT4l) . It is reasonable then to compare values of Eog in simulations SI -S3 with those of E,„„v obtained in other 
experimental and numerical settings. In particular, in S3, E0.5 = 0.73 mm"'. This is comparable to the range of values 
of Emflt w 0. 12 - 0.6 mm"' obtained in a number of experimental studies using a wide variety of flame configurations 
(see 1I3I for a summary table). The fact that the value in S3 is somewhat larger is not surprising, given that the turbulent 
intensity in it is substantially higher than in the experiments. 

Overall, however, the traditional definition of E [3] is not equivalent to the definition of Ey used here unless 
F = 0.5. Therefore, unlike E„„v, E0.5 is not necessarily the maximum value of Ey for all values of F. Fig. |4f shows 
that while in SI E05 is indeed approximately the largest value, this is not the case in S2 and S3. Nonetheless, in all 
three calculations, E0.5 is representative of the values of Ey throughout the reaction zone with the difference between 
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!E().5 and S0.15 being ^6-9% in the higher-resolution cases. 

3.3. Relation between Ay and 

Fig. |5] shows the relation between the isosurface area and density in three key regions of the flame: the peak 
reaction rate {Y - 0.15), the boundary of the preheat and reaction zones (Y - 0.5), and the coldest region of the 
preheat zone (Y = 0.95). A pronounced linear correlation between Ay/L^ and Ey exists throughout the flame interior. 
This demonstrates that at all values of Y, an increase of Ay is greatly influenced by the associated tighter packing of 
the isosurfaces inside the flame brush. 

With increasing Y the average Ay 1 1? and Yy become larger, as was seen in Fig.|4^,f, due to the overall shift of 
the distributions toward larger values of both Ayjl} and Ey. Note also a slight increase both in the scatter of the 
distribution and the slope of its least squares fit at higher Y. Thus, isosurfaces of higher fuel-mass fractions are indeed 
consistently more tightly packed than at lower Y. This is in agreement with the observation that the flame surface is 



wrinkled on finer scales on the fuel side than on the product side 111 111 . 



3.4. Distributions of Ay and Zy and the effects of small-scale turbulence 

How do these results concerning the distributions of Ay and Sy, shown in Fig.|4^,f, relate to previous conclusions 
ifTT l regarding the inability of the turbulent cascade to penetrate the flame interior and the role of small-scale turbulent 
motions? As was discussed in § l2.1l (also see 1, 1 11I for further details), the progressive increase in resolution from SI 
to S3 extends the turbulent cascade to smaller scales, and this leads to a substantial increase in the energy of turbulent 
motions on scales A < Si in nonreactive turbulence. This increase in resolution was found to cause the flame surface 
on the fuel side to be wrinkled on progressively finer scales, while remaining virtually unchanged on the product 
side. Furthermore, it was determined that the deviation of the internal turbulent flame structure from that of the 
planar laminar flame increases with decreasing temperature. These two results showed that the eff'ects of small-scale 
turbulent motions are most pronounced in the coldest parts of the preheat zone. With increasing temperature, these 
effects diminish and they completely disappear once the reaction rate becomes significant [11]. 

The results presented above are consistent with this picture. The steady growth of Ay/L^ and Ey on the fuel 
side of the preheat zone with increasing resolution shows that more intense small-scale motions are indeed present 
there, folding isosurfaces on progressively finer scales. It must be emphasized that the change in Ay alone does not 
necessarily mean that it is caused by finer wrinkling on smaller scales. Only when viewed in conjunction with the 
surface density, which shows a very similar dependence on Y, does such increase in Ay serve as a strong indication of 
the small-scale wrinkling. 

With increasing temperature, not only do Ayjl} and Ey decrease rapidly, but also profiles for S2 and S3 begin 
to approach each other This means that the main difference between these two simulations, namely the presence of 
more energetic small-scale turbulence in S3, is being eUminated and, thereby, motions on scales A < 6l are gradually 
suppressed. Since it is these small scales that enhance the diffusive transport, which in turn broadens the preheat zone. 



their suppression causes the internal flame structure to approach that of the planar laminar flame 111 111 . 

As the reaction rate becomes substantial at F < 0.5, both Ay/L^ and Ey become similar in S2 and S3. This 
suggests that, at this point, only scales A > 6l remain energetic. Indeed, these scales are originally the same in both 
calculations (cf. Fig.[T]). Consequently, they generate similar isosurface areas and densities. Since these scales cannot 
supp ort small-scale diffusive transport, any broadening of the internal flame structure effectively disappears at 7 < 0.5 
11 111 . Further decrease in Ay/L^ and Ey is the consequence of the continued suppression of the progressively larger 
scales A > 6l. 

There is one important distinction between the distributions of Ay/L^ and Ey and the time-averaged profiles of Y 
and T, which represent the internal flame structure fTT). The Y and T profiles showed very little variation between 
simulations SI -S3 not only in the reaction zone, but also in the preheat zone. This is in contrast with the distributions 
Ayjl} and Ey, which differ substantially in the preheat zone even between S2 and S3. Such discrepancy suggests that 
small scales, A < 6l, contribute differently to the turbulent diffusive transport and to the wrinkling of the isosurfaces. 
The former is primarily governed by the largest of these small scales, i.e., scales not much smaller than 6l, since these 
scales are associated with the highest velocities. The energy contained on these scales is close in all three simulations 
and, thus, they produce the same structure of the broadened preheat zone. At the same time, all scales smaller than 
6l contribute to the isosurface wrinkling, and the different energy content of these scales causes the distributions of 
Ayjl} and Ey to differ substantially in the preheat zone. 



9 



Qualitatively, both the internal flame structure described in Ol ill and the distributions of Ay and presented here 
create a consistent picture of the transformation that turbulence undergoes as it passes through the flame. Changes 
in the energy budget between different scales, which are most certainly accompanied by the development of both 
anisotropy and inhomogeneity of the velocity field, are complex. While the evidence presented here and in 111 1 1 
provide hints about the nature of this transformation, the details remain unknown. For instance, are the energy release 
and the resulting fluid expansion the only effects responsible for altering the turbulent field and redistributing energy 
between different scales? What are the relative contributions of various small scales to the increase in Ay and Zy? 
How does the shift in balance between small and large scales with increasing temperature change at different turbulent 
intensities? All of these questions need to be addressed in future studies. 



4. What is the flame surface area? 

In the wrinkled and corrugated flamelet regimes, in which eq. ([TJ is typically applicable, the flow is laminar on 
scales A < 6l. As a result, isosurfaces of different values of Y are parallel to each other, the flame is folded by 
turbulence as a coherent structure, and, therefore, it has a well-defined surface area At- Consequently, the stretch 
factor / has a unique value at each moment in time, and so it directly shows how much of the increase of St can be 
ascribed to the increase of A 7-. 

In contrast, in the thin reaction zone regime discussed here, results presented in § [3] demonstrate that the flame, 
folded inside the flame brush by high-intensity turbulence, has a complex internal structure. The various parts of the 
flame have distinctly different responses of its various parts to the action of turbulence. Consider, by analogy with 
eq. dU, the factor ly, which relates the increase of the turbulent flame speed relative to its planar laminar value to the 
increase of Ay rather than At, 

In this case, unlike /, ly is a function of Y. In particular, based on the data shown in Fig.|4] ly can vary by as much as 
a factor of 4 depending on the choice of Y. As a result, values of ly close to and substantially larger than unity can be 
found in different regions of the flame. 

This shows that before we can answer the question whether the observed S t can be accounted for by the increase 
in the flame surface area, we must first determine what this area is in such high-speed turbulent flows, and whether 
such a concept is even applicable in this regime. In particular, what value of Y fully and accurately represents the 
overall behavior of the turbulent flame? We address these questions by considering the correlation of Ay with two 
global properties of the turbulent flame, its speed and width. 

4.1. Relation between Ay and S t 

Fig.|6]shows correlations between St IS l and Ay/L^ calculated for the same three representative values of Y as in 
Fig E] namely Y - 0.15, corresponding to the region of peak reaction rate, Y - 0.5, corresponding to the boundary 
between the reaction and the preheat zones, and Y = 0.95, corresponding to the coldest part of the preheat zone. Left 
panels of Fig. |6]show the time-evolution of all quantities, and the right panels show the corresponding correlation 
scatter plots. 

Since both the local flame speed and the induction time of the unburned fuel have finite values, there must exist 
a delay in the response of St to the changes in the flame configuration, represented by the areas of its isosurfaces. 
Therefore, the degree of correlation between St and Ay is a function of a time lag, Af, between these two quantities. 
For each value of Y in Fig. |6] we determined the time lag, Af^. , that produced the best correlation. This was done 
by directly calculating the cross-correlation between ST{t)IS i and Ay(t)/L^ and finding its maximum. The value 
obtained for At^ was verified by determining the least squares fit for the distribution of STit)/S l vs. Ay(f - At')/L^ 
and by finding At'^ that maximized the slope of the fit and minimized its residuals. The values of Ay/L^ in Fig.|6]were 
then shifted in time with respect to St/S l by the corresponding time lag Ate given in the lower right corner of the 
panels (a), (c), and (e). 

The first key conclusion emerging from Fig. |6]is that St and Ay become progressively less correlated with in- 
creasing Y. While the correlation is pronounced throughout the reaction zone and is exceptionally strong near the 
peak reaction rate, S t and Ay are only very weakly correlated in the colder parts of the preheat zone. Note that panels 
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(b), (d), and (f) have the same scale, which shows the relative increase in the scatter of the distribution and its overall 
shift to higher values of Ay in the preheat zone. 

The time delays that give the best correlation between the two quantities are positive, and thus they show that S t 
does indeed lag behind Ay- Moreover, the magnitude of Ate increases with increasing Y. Qualitatively, this agrees 
with the fact that St 31 any given moment is primarily determined by the region of the highest reaction rate, which 
corresponds to the values of Y close to 0.15. Consequently, the time delay between St and the area of the isosurface 
representing the peak reaction rate is close to zero, namely O.QATed- Since turbulence reorganizes the flame on a 
timescale ~ Tgd, such a small lag is too short for the flame structure to change in any significant way. Therefore, the 
correlation between Ao. 15 and 5r is very strong. 

At the same time, since the reaction rate is fairly low at F = 0.5, the contribution of this region to the overall 
turbulent flame speed is small. Burning still needs to accelerate substantially in this area in order for it to become the 
new region of peak reaction rate, which will then predominantly determine the magnitude of St- The time required 
for this to occur results in a longer time delay, namely QAlzed- During this time, however, turbulence is able to change 
the flame structure more strongly than in the case of F = 0.15. As a result, the degree of correlation between Aq 5 and 
St decreases. 

In the coldest regions of the preheat zone, there is no correlation between S t and Aq 95 at small values of At^. We 
found evidence of weak correlation for a much larger time delay At^ - 0.6 It^,/. This value, however, is comparable 
to the time during which the turbulence completely reorganizes the structure of the turbulent flame. Therefore, by 
the time burning reaches the reactants in the outer regions of the preheat zone, the connection between their original 
distribution and the resulting turbulent flame speed is significantly disrupted. 

We also found weak correlation between A0.95 and A0.15 with the time lag for A0.95 being Ate = 0.56Ted and the 
slope of the least squares fit of 0.49070 This weak correlation between A0.15 and A0.95 is consistent with a similarly 
weak correlation between S t and A0.95 found at a comparable time lag. In particular, it reflects the fact that after the 
time x;0.6Ted, burning reaches the outer part of the preheat zone. As a result, the former Y - 0.95 region becomes the 
new Y - 0.15 region and thus the new site of the peak reaction rate which now determines St- Substantial change in 
the flame-brush structure during this time, however, again leads only to weak correlation between Aq 15 and Aq 95. 

In order to verify this physical meaning of the obtained time delays, values of Ate can be compared with the 
characteristic induction time, T,„rf, in the planar laminar flame. The T,„rf is the time for a fluid element to reach the 
peak reaction rate starting at some initial temperature or, equivalently, fuel mass fraction. In reactive mixtures with 
large activation energies, such as Ha-air considered here, t,„j is also approximately the time necessary to reach the 
post-flame temperature, Tp ll36il . 

The Tind can be calculated using the velocity distribution Ul{x) in the exact planar laminar flame solution as 

Ti„d(Y')^ (17) 

Jx{Y') Ul{x) - f/oo 

Here Uoo is the fuel velocity at infinity, Y' is the initial fuel mass fraction of the fluid element, and Y* is the fuel mass 
fraction corresponding to the peak reaction rate. Instead, to maintain consistency with the numerical simulations, we 
computed Ti„d directly by solving the full time-dependent reactive Navier-Stokes equations O)-® on a uniform ID 
grid with the same numerical integration method used for S1-S3 (§ 12. 11 1. The flame was initialized as a discontinuity 
with the fuel and product temperature and density equal to To, po and Tp(), ppfi, respectively, given in Table [S] 
Constant pressure Pq and zero velocity were set throughout the domain. We introduced an additional advection 
equation describing a scalar ^ that was not subject to diffusion or reactions and was only advected with the flow. After 
the steady-state structure of the flame had been established in the computational domain, ^ was set to unity in one cell 
corresponding to the desired initial T' or, equivalently, Y', and to zero everywhere else. As the system evolved, the 
fluid element marked with the scalar ^ moved through the flame and burned. The time it took this fluid element to 
reach peak reaction rate, i.e., Y - 0.15, gave r,„rf(y')lj 



'in contrast, we were unable to find any statistically significant correlation or anticorrelation between So. 15 snd Z0.95 at any time lag. 

'^The distribution of f would spread with time over several cells due to numerical diffusion. To compensate for this, we used very high resolution, 
namely Ax = Slq/HS, which is much larger than the resolution necessary to obtain an accurate laminar flame solution, i.e., typically Ax = Stf)!^ 
t32ll . As a result, the distribution of f was found to spread over the range of values of Y with the full width at half maximum equal to AF a 0.05. We 
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Using this procedure, we found that for a planar laminar flame, 



T,w * 0.27t„/ * 2.46Af„ for 7 = 0.5, 

T,„rf ^ 1.29t„/ ^ 2.12Afr, for Y = 0.95. ' 



When turbulent heating of the fuel is taken into account (see § 12.2b . i.e., when the fuel temperature in this calculation 
is increased to 320 K from Tq = 293 K while maintaining constant density po, these values of r,„^/ decrease by 
~ 20 - 25%. Further details and, in particular, the full distributions of r,„^/ throughout the flame for these two fuel 



temperatures, can be found in Appendix B 



The calculated time lags are, therefore, within a factor of ^2 of the corresponding values of T,„rf, which, given the 
statistical nature of Afc, is reasonable agreement. This shows that the correlation time delays for the distributions of 
St/S L and Ay/L^ in the turbulent flame can indeed be associated with the induction times necessary for the burning 
to accelerate in regions with higher reactant concentrations. 

The values of Ate, however, are substantially lower than the characteristic propagation time of the laminar flame, 

Slo 30 L 

Td = — = = 3J5Ted, (19) 

Sl,o 8 t/ 

where we used the values of L and U given in Table |3] Even for Y = 0.95, At^ = Q.6lTed ~ (1/6)tj. Furthermore, 
this Ate is an order of magnitude smaller than the time l-67SLfl/S l,o ~ 6.3T(.rf, which is necessary for the flame to 
propagate over the characteristic distance 1.67(5l,o separating Y = 0.95 and Y - 0.15 in the laminar flame structure 



(cf. Fig. 7 in nun ). 

A good analytical estimate for Afc throughout the reaction zone can be obtained using the expression for the 
adiabatic induction time ll36ll . 

where T' and p' for a given Y' are determined based on the exact planar laminar flame solution, and the values of B, 
q, and Q are given in Table |2l In particular, at the fuel temperature To and density po, 

t«^^0.05t,^* 1.25Af„ for 7 = 0.15, 

-rf^^ * O.llTed ~ 1.91Afr, for Y = 0.5. ^ ' 

Using the specific heat at constant volume C,. - R/M(y - 1) instead of Cp in eq. (l20t . i.e., assuming that burning 
occurs at constant volume rather than at constant pressure, would decrease the values of r"^^ by w 15% to O.QATed 
and 0.1 8Te^, respectively. This shows that, even though is typically applicable in the context of an autoignition 
process rather than propagation of a flame, for Y - 0.5 it provides a somewhat better approximation to the calculated 
value of Ate than t,„£/ given in eq. (fTSl l. Further details, including the effects of turbulent fuel heating on t^!^^, are given 
inf " ^ 



Appendix B 



Finally, in addition to the strong correlation between the values of St and A0.15, the distribution shown in Fig.|6j3 
has another important property that is not present at higher values of Y. The linear least squares fit for the case 
Y = 0.15, shown with a solid line in Fig.|6j3, has the form 

|l = 1.17^ _ 0.08. (22) 
Sl 

For Ao.is/L^ = 1, this expression gives St/S l - 1.09. This demonstrates that as the flame becomes less convolved, 
i.e., as AoAs/L^ — » 1, it behaves progressively more like a planar laminar flame which is manifested in St — * St. 
Linear least squares fits for the two other values of Y, shown in Fig.|6}l,f, do not recover the value of St/S l = 1 as 
Ay/L^ 1. 



assumed that the position of the fluid element being tracked was marked by the maximum of this distribution. A convergence study with increased 
resolution, which caused f to be distributed over progressively nan'ower range of Y, was used to verify that the obtained values of t,„^ were indeed 
accurate. 
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4.2. Relation between Ay and 6t 

Next consider the correlation of A 7 with the second key global characteristic of the turbulent flame, namely its 
width, ^T". Fig. |7] shows the correlation between Ay/L^ and 6t/L for the same three values of Y as in Fig. [60 The 
distribution of 6t/L as a function of AyjL^ shows the same trend observed for St/S l- In particular, the correlation 
between the two quantities is also the strongest for Y - 0.15 and decreases with increasing Y. The values of 6t/L 
and A0.95 appear to be almost completely uncorrected. Note, however, that the correlation in the reaction zone, and, 
in particular, in the region of peak reaction rate, is weaker for 6t than St, with the distributions of 6t having a much 
larger scatter 

In order to illustrate the relation between the three key quantities characterizing the instantaneous flame-brush 
structure, namely Ay, dr, and I.y, each data point in Fig.|7]is colored according to the corresponding value of "Ly (cf. 
Fig. |5]l. At lower Y, the increase of the isosurface area is associated both with the increase in the overall extent of 
the flame brush and with tighter packing of the isosurfaces. With increasing Y, however, the correlation of Ay/L? 
with 6t/L decreases while remaining quite pronounced with Hy. This shows that the area of the isosurfaces of higher 
Y changes primarily only due to their wrinkling on progressively smaller scales, as discussed in §[3] (also see jl in ). 
Indeed, Sj is determined predominantly by the flame folding on the largest scales. Consequently, the weak correlation 
between Ay and 6t at high Y demonstrates that the increase of the areas of these isosurfaces cannot be associated with 
such large-scale wrinkling. Instead, this increase must be translated primarily into tighter packing of the isosurfaces 
causing Ay to be correlated with Y.y. 

5. What determines St in the thin reaction zone regime? 

5.1. Can St be explained by the increase in the flame surface area ? 

The strong correlation of A0.15 with S t and 6t, which represent both the global energetics and the global structure 
of the turbulent flame, show that, in the thin reaction zone regime, isosurfaces of Y close to the peak reaction rate 
accurately characterize the overall evolution of the turbulent flame. Consequently, for the reactive mixture considered 
here, the flame surface area. At, and density, E7, are represented by the isosurface ofY- 0.15. 

This conclusion is further supported by the limiting behavior of the distribution of St/S i, which approaches unity 
as Ao.is/L^ — > 1. In contrast, atY- 0.5, St/S l becomes progressively less than Aq.s/L^ at smaller isosurface areas. 
Furthermore, at F = 0.95, almost afl values of St/S i lie below the St/S l = Ao.gs/L^ line, and they become almost 
half of Ao.95/L^ at larger isosurface areas. In Le = 1 reactive mixtures, there are no mechanisms that can lower the 
local flame speed below S l, which would be necessary to produce values of St/S l lower than the increase of the 
flame surface area. Consequently, such properties of the distributions of St/S l confirm that isosurfaces of higher Y 
cannot represent the flame surface. 

With this definition of the surface area of the turbulent flame, we can now revisit the question of the relation 
between At and the speed of the turbulent flame in the thin reaction zone regime. This relation can be represented 
by eq. dTSI ) in which Y — 0.15, similarly to eq. ([TJ applicable in the wrinkled and corrugated flamelet regimes^ 
Fig.[8]shows the ratio /0.15 - {St/S l)/{Aq,\s/L?') for the distribution given in Fig.|6j3. Time-averaged values /0.15 - 
(S T / S i) / (Aon / L?-) for all three calculations are listed in Table I4f1 In particular, 7o.i5 = 1.14 in S3, and this value can 
be viewed as converged to within a few percent with a faster than linear rate of convergence. 

These results lead to the following two key conclusions: 

1. In the thin reaction zone regime, the instantaneous turbulent flame speed is also primarily determined by the 
increase of the flame surface area. 



'since both Ay and Sr represent the same instantaneous configuration of the flame, no time delay between them is expected. Consequently, no 
time shift was applied to either quantity in Fig.0 As a check, however, we verified that indeed the cross-correlation had its maximum at the zero 
time lag. 

^In light of the fact that 5 7- and Ay are best correlated when Ay is shifted in time with respect to 5 7-, the question arises whether such a time 
shift should also be applied in eq. )161 . In principle, an argument could be made that, due to the inherent delay in the response of 5 7- to the changes 
in the flame structure, the current value of 5 7- is the result of Ay that existed earlier in time and, thus, this earher value must be used to assess the 
relation between 57- and Ay. At the same time, as was shown, the time lag At^ is extremely small for Y = 0.15. During this time, the change in the 
value of Ay is negligible and applying the time shift in eq. )16t does not cause any appreciable change in the result. 

'Using instead the expression I0.15 = (S r / S l)/{Ao.is / L~) gives virtually the same result. 
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2. In addition, in the course of system evolution, S t periodically exhibits an exaggerated response to the increase 
in At- For the system size and turbulence intensity considered here, this excess burning velocity can become as 
high as 30% of Ar/L^. 

This exaggerated response in simulation S3 is represented by values of /0.15 > 1 in Fig. [8] and is illustrated in 
Fig.|6^,b as the shaded gray area. Such a substantial deviation of /0.15 from unity brings up the following question: 
what mechanism is responsible for raising S t beyond what can be attributed to the increase in A7? 

5.2. Potential causes of the excess values of St 

As was discussed in § [1] Damkohler's concept [T] implies that the burning velocity of a turbulent flame in the 
presence of sufficiently intense turbulence is determined by two processes: wrinkling of the flame surface by large- 
scale turbulent motions increasing At, and the enhancement of the local flame speed S„ by the small-scale turbulence 
providing stronger diffusive transport. While the simulations show that the effect of large-scale turbulence is indeed 
present and has a dominant effect on St, can the observed excess values of 5 r be ascribed to the action of small-scale 
motions? 

Consider the following important characteristic exhibited by the distribution of /0.15. Along with /q 15 substantially 
larger than one, the system also periodically develops values of /q 15 within only a few percent of unity, which is 
comparable to the uncertainty in estimating the instantaneous Aq 15 and, thus, /0.15 (see also Fig.|6h). This shows that 
at those instances when /0.15 ~ 1, S,, at all points of the flame surface must be almost identical to S l. Consequently, 
the small-scale turbulence entering the flame, if it were directly determining S„, would have to periodically undergo 
transitions between almost complete suppression, when /0.15 ~ 1 and S„ ~ S l, and significant enhancement, when 
/0.15 becomes as high as 1.3. 

While we are unaware of any potential mechanisms that could cause such periodic complete suppression of small- 
scale turbulence in the cold fuel, our analysis also does not show any changes in the state of turbulence that could 
serve as a manifestation of such suppression. In order to determine whether there exists a connection between the 
instantaneous turbulent field in the domain and the resulting burning speed, we have analyzed the correlation of 
several turbulence characteristics with both St and /0.15. In particular, we considered the total rm.s. velocity in the 
domain. Furthermore, to characterize turbulence inside the flame brush and in the fuel entering it, we also considered 
individual components of velocity, m,, as well as kinetic energy, 0.5pu^, averaged over the volume of the flame brush 
and the cubic region of size L located in the fuel immediately ahead of it. We were not able to find virtually any 
correlation between any of these quantities and either St or /0.15. As an example. Fig. |9] shows the representative 
correlation scatter plot of /0.15 as a function of the normalized x-component of kinetic energy Ek,x = (pu^)/{poUf) 
volume-averaged over the region discussed above. Similarly to the case of the distribution of /0.15 vs. Aq.is/L^ in 
Fig.|8] the increasing trend seen in the least squares fit is too weak to determine its statistical significance based on the 
available data. 

These considerations suggest that the enhancement of S„ over S ihy turbulent transport does not provide a plau- 
sible explanation of the excess values of 5 r periodically developed by the flow. Consequently, a different mechanism 
must augment the effect of the increased At- 

Instead of being the result of a uniform global increase of S„ throughout the flame brush, accelerated burning can 
be caused by a significant increase in 5 „ in isolated regions of the flame surface. In those regions, however, S„ must 
be much larger than S lin order to significantly raise the overall burning speed of the turbulent flame. 

Such local enhancement of S„ cannot be caused by the local inhomogeneities in the thermodynamic state of the 
fuel0 In particular, in a subsonic turbulent flow, such as the one considered here, with the turbulent Mach number 
Maf substantially less than unity (see TableO, local variations of density or temperature are much too small to affect 
S„ in any appreciable way. For instance, it requires an increase of fuel temperature of over 150 K to double S„. We 
do not observe temperature fluctuations of such magnitude in the flow. 

S n can also be increased by the substantially higher local turbulence intensity associated with intermittency. While 
we do observe intermittency in the flow field in our simulations, regions of large velocity enhancement are statistically 
too rare, their spatial extent is too small, and they are much too short lived to increase S„ and, thus, to make any 



**The increase of S„ due to the global fuel heating by turbulence was taken into account in St, as discussed in § 12.21 
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significant contribution to the overall increase in S't-. Potential effects of turbulence intermittency on the flame have 
been studied by Pan et al. fjT] in the context of the Rayleigh-Taylor-driven thermonuclear flames in the interior of 
a white dwarf during a Type la supernova explosion. Their results suggested the possibility that intermittency in the 
velocity field could locally disrupt and broaden the flame, which potentially would increase the local burning speed. It 
is not immediately clear, however, to what degree their conclusions apply to the system discussed here. The conditions 
in the white dwarf interior considered in 137,1 are substantially different from those in S1-S3. Moreover, their analysis 
did not consider any feedback of the flame on the turbulent field and, in particular, it did not include fluid expansion 
due to heat release which invariably affects intermittency of the flow. At the same time, the role of intermittency in 
the dynamics of turbulent flames does require further investigation. In particular, it would be important to extend 
the analysis of Pan et al. OTIl by including the feedback of the flame on turbulence under conditions characteristic of 
chemical rather than thermonuclear flames. 

Instead, we suggest that significant local increase of S „ does indeed take place in the turbulent flame due to the 
flame collisions resulting in the creation of regions of large flame curvature. We discuss this process next. 



6. Accelerated burning caused by flame coUisions 

6.1. Local increase of S „ in cusps 

According to the theory of flame stretch (see ff\ for a review), when Le - I, the flame is not affected significantly 
llssll either by flow-induced strain or curvature. In particular, a stationary spherical flame supported by a point source 
of mass, which is an example of a curved unstrained flame, has an internal structure, and thus local burning velocity, 
identical to that of a planar laminar flame ff\. Such analysis, however, is typically carried out for a flame with the 
curvature radius rc » Sl- Using the physical model and the Athena-RFX code described in § 12.11 we also performed 
simulations of an idealized spherical and cylindrical flame propagating inward into the stationary fuel. The results 
of these simulations are, in fact, consistent with the theory. Furthermore, they show that S „ - S l not only when 
fc » ^L, but up until the moment when the flame curvature becomes xI/6l, i-e., until the flame effectively collapses 
onto itself. At this point, S„ increases substantially. Beyond this moment, however, the flame itself effectively ceases 
to exist and, therefore, the local flame speed looses its meaning. Such an idealized situation is not representative of 
the conditions arising in an actual turbulent flame. Nonetheless, regions of large flame curvature ~ 1 /5l are frequently 
created in the high-speed turbulent flow considered here, and they naturally provide a mechanism for a significant 
local enhancement of S„. 

The results presented in §|3]showed that intense turbulence causes tight folding and packing of the flame inside the 
flame brush. Consider the inverse of the surface density, which is a measure of an average separation between surface 
elements. Given that E0.15 = 0.67 mm"' in simulation S3, the average distance between individual flame sheets is 
l/So.is ~ 1-49 mm « A.lSi. At the same time, Eq is can be as high as 0.97 mm"' (Fig.|5h), resulting in an even lower 
separation of ki 3>.2Sl- This is comparable to the full flame width If ^ 26l (cf. Fig. 7 in 1 11]). Such tightly folded 
flame configurations invariably result in frequent collisions of individual flame sheets. 

Fig. [To] gives an example of such collision in S3. The figure shows the flame-brush structure at three times: 
11.86Tf^ (upper panel), as well as Q.lTed = 3 fis and Q-lxed = 6 //s later (middle and lower panels, respectively). 
Initially, a region with a highly convolved flame develops in the flame brush (region A). At this time, S„ ~ S l locally 
throughout region A and changes in the flame configuration due to self-propagation can be neglected on timescales 
considered in Fig.[TO] As turbulent motions continue to bring individual flame sheets closer to each other, the curvature 
radius of the flame becomes close to 6l, and the preheat zones begin to overlap substantially over an extended region 
of the flame surface. This marks the formation of two regions of high flame curvature ~ 1/6l (regions B), which 
we will refer to as "cusps." Regions C in the lower panel show that 3 fis later, the flame sheets have merged and 
formed two extended reaction zones, which suggests substantially accelerated burning in that area of the flame brush. 
Note also a rapid decrease of the area of the Y - 0.5 isosurface represented with the thin black line. Comparison of 
regions B and C shows that this isosurface propagated over the distance a 0.5 mm over the time 3 yus, which implies 
the propagation speed (but not the local burning speed) ^ 555 l. During the time shown in Fig. [10] two other highly 
elongated regions of flame collision have formed near the upper face of the domain. Thus, Fig.[TO]illustrates that flame 
collisions are ubiquitous in the turbulent flame, and they do produce regions of large flame curvature. Moreover, their 
evolution in Fig.[Tn]suggests substantial increase of both the local burning and propagation speeds in the cusps. 



15 



Qualitatively, two different channels of cusp formation can be identified. They are shown schematically in Fig.fTTI 
In the first case, turbulent motions on a certain scale A with velocity Ua stretch and fold the flame. As a result, a 
narrow elongated structure is formed. When two flame sheets effectively collide, the curvature of the flame surface 
becomes ~ 1/6l and the cusp is formed with a characteristic length Ic ~ A. We will refer to this channel as "local 
flame collisions," since the cusp is formed by flame regions located close to each other both in 3D space and along the 
flame surface. In the second case, turbulent motions fold the flame on a larger scale, thus not increasing its curvature 
substantially. Two cusps are formed when two smoothly curved flame sheets collide. We will refer to this channel 
as "nonlocal flame collisions," since the colliding flame regions can be located far from each other along the flame 
surface and they can, therefore, originate from different regions of the flame brush. In this case, the characteristic 
length of the resulting cusp, l^, can be significantly smaller than the scale. A, of the motion that produced the cusp. 

In an actual turbulent flame, cusps are often created under the combined action of these two channels, an example 
of which is seen in Fig. [10] Region A is initially formed by the flame surface being stretched and folded on a large 
scale, as in Fig. [TTb . At the same time, smaller scale motions fold the flame surface locally, as in Fig. [TTh . forming 
smaller cusps seen in regions B. 

While both local and nonlocal flame collisions are equivalent in terms of the resulting cusp properties, they differ 
in one important respect. In the first case, curvature of the flame surface at a specific point increases gradually until 
it becomes large marking the formation of a cusp. In the second case, the curvature remains small at all points of 
the flame surface until the actual moment of collision when two flame sheets merge and the region of large flame 
curvature is abruptly formed. These distinctive characteristics of the two processes will be important for the analysis 
of their relative efficiency to form cusps at a given turbulent intensity and system size. 



6.2. Structure and properties of cusps 

Formation of cusps (or "folds") by the wrinkled flame in a turbulent flow was first suggested by Karlovitz et 
al.lH^ (also see discussion in f4d\). Their dynamical role was first considered by Zel'dovich [41] (also see [i4Zl ). who 
suggested them as a mechanism of flame stabilization that prevents the unbounded exponential growth of the flame 
surface under the action of the Landau-Darrieus instability |43]. Such regions of large curvature would have a high 
propagation velocity, causing rapid annealing of the flame surface, and thus would provide an efficient mechanism to 
counterbalance flame instability |41]. It must be noted that both Karlovitz et al. ll39ll and Zel'dovich |41] assumed that 
cusps form due to the self-propagation of a wrinkled flame surface following the Huygens principle, rather than as a 
result of flame collisions by turbulence. 

It was suggested early on both by Lewis and von Elbe||40ll and Zel'dovich 04111 that, due to the focusing of the heat 
flux into fuel, S„ ceases to be equal to S l near th e tip of the cusp, and the cusp region develops a structure similar to 
the tip of the Bunsen flame. The flame in |39, 4^ 41 , i^l, however, was considered to be a gasdynamic discontinuity, 
and so its internal structure was ignored. As a result, the actual increase of the local flame speed in the cusp could not 
be determined. 

In order to analyze quantitatively the effect of flame collisions and to determine the actual increase of S„ in the 
resulting cusps, consider an idealized model of the cusps observed in Fig. [10] In particular, consider two symmetrically 
located planar flame sheets propagating with the speed S la and approaching each other with the inclination angle to 
the mid-plane a. At the point of collision, the flame sheets merge and form a cusp. Its properties, such as its structure, 
speed of propagation, and effective burning velocity, are determined in this configuration only by two parameters, 5l,o 
and a. Such cusps can be formed through either of the channels shown in Fig.fTTI 

Fig.[T2]illustrates the flame structure formed in this situation. It shows the distribution of the fuel mass fraction in 
two simulations performed for a = 1° and or = 4° with Athena-RFX using the same physical model as in simulations 
S1-S3. The domain has the resolution Ax - <5l/32 and zero-order extrapolation boundary conditions on all sides. At 
f = 0, two intersecting flame sheets were initialized with the exact structure of the planar laminar flame corresponding 
to the fuel temperature Tq and density po (see Table|2]i. Uniform pressure Pq and zero velocities were set in the domain 
att - 0. After the initial transient stage, the flame develops the structure shown in Fig. [12] The structure of the cusp 
itself does not change with time, provided that the planar flame sheets extend sufficiently far from it. Therefore, 
Fig. [T2]can be viewed as a steady-state solution. 

Fig. [12] shows that flame collision at very low inclination angles results in the formation of a highly elongated 
structure. It is formed by rapid heating and subsequent ignition of a larger amount of fuel than in the planar laminar 
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flame due to the focusing of the thermal flux from two approaching flame sheets. Consequently, the extent of this 
structure is determined by the size of the region in which the preheat zones of approaching flames overlap and, thus, 
create the necessary enhanced thermal flux. Since the width of the preheat zone is ^ 6^^, very small values of a are 
required for the cusp to be substantially broader than 6l^[). This can further be seen in the distributions of Y, T, and 
Y along the symmetry axis of the cusp shown in Fig. [13] The structure of the reaction zone in the cusp very quickly 
approaches that of the laminar flame with increasing a, so that akeady at a ~ 4° the two reaction zones become very 
similar. At the same time, the preheat zone of the cusp remains substantially broadened for much larger values of a. 

There are three characteristic speeds in this problem. Sufficiently far from the tip of the cusp, the flame locally 
propagates normal to its surface with the laminar flame speed 5l,o, i-C-, there S„ = 5l,o- The cusp itself moves into 
the fuel with the phase velocity Uc, which results simply from the two inclined planar surfaces colliding with each 



other This speed can be determined based on purely geometrical considerations 114111 . 



= (23) 
sma 

This cusp propagation speed, normalized by S l,o, is shown as a solid line in Fig.[T4h along with Uc/S u) determined 
as the velocity of the leftmost point of the F = 0.15 isosurface in simulations for four values of a. Eq. ( l23T l is within 
< 1% accuracy of the computed values. Uc approaches the laminar flame speed as a 90°, i.e., when the two flame 
sheets cease to advance toward each other. At the opposite limit of small a, Uc can be quite large and, in principle, 
it could be infinite when a ^ 0°. In particular, at o- 0.5°, Uc becomes larger than the sound speed in cold fuel 
indicated with the horizontal dashed line in Fig.[T4h. Note also that at o- = 1°, the value of Uc ~ 57S la is comparable 
to the high propagation velocity «555l of cusps in regions B - C in Fig.fTOlestimated in § 16.11 

The third characteristic speed, which is the most important for us, is the local burning speed in the cusp, S„. The 
broadened reaction zone at low flame-inclination angles (Figs. [T2]and[T3b) causes more fuel to be consumed per unit 
flame surface area than in the planar laminar flame, and, therefore, it leads to a larger local flame speed S„ > 5l,o- 
Since the reaction-zone width is the largest at the tip of the cusp, S„ has its maximum there and it gradually decreases 
to its laminar value 5l,o in the planar regions of the flame sufficiently far from the cusp. The maximum value of S„ at 
the cusp tip can be found as 

S: = max(5„) = - f pYdx = -- f p^Fexpf - ^] dx. (24) 

Here eq. ^ was used for Y, and the integral is taken along the symmetry axis of the cusp shown in Fig. [12] Fig.[T4b 
shows the computed values of 5* normalized by S 1,0 for the same four inclination angles as in Fig.[T4h. 

Two key conclusions emerge from Fig. [14] First, 5* is substantially lower than Uc, being less by more than an 
order of magnitude at the values of a considered here. Second and most important, 5 * can indeed be much larger than 
Sla at small inclination angles. 

6.3. Connection between the increase of S „ in cusps and St 

We are interested, however, not just in S„ but, rather, in the burning speed of an extended flame region, as this is a 
direct equivalent of the turbulent flame speed. The total burning rate of the flame configuration shown in Fig.[T2lis 

niR 

Sc = — , (25) 

PO 

where rriR is the total mass of reactants converted into product per unit time. If the flame were to propagate everywhere 
with the speed Sla, then Sc - AyS u), where Ay is the surface area of the flame in such a configuration. In this case 
all isosurfaces of Y would be paraUel to each other and, thus. Ay based on any value of Y could be used. Otherwise, 

Iy^-^>1. (26) 

AyC) L.O 

Here ly is completely equivalent to the definition used in eq. dTSI l in the context of the turbulent flame speed. Indeed, 
substituting the definition of St given by eq. (ITTT l into the definition of ly given by eq. (ITSI l shows that mrIpo — 
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hAyS L.o = Sc- Furthermore, as can be seen in Fig. [12] at low a isosurfaces of Y do not have the same surface area. 
Therefore, similarly to the case of the turbulent flame, here ly is also the function of Y. For consistency with previous 
results, we will consider /0.15 based on the isosurface of F = 0.15 representing the peak reaction rate. 

Since S„ monotonically decreases away from the cusp, both S c and /0.15 depend on the size of the flame region 
being considered. This size can be represented by the variable length 1^, illustrated in Fig. [12] which is defined as 
the distance from the leftmost point of the isosurface of F = 0.15 in the direction of cusp propagation. Considering 
different values of Ic is equivalent to varying the fraction of the total flame surface area in Fig. [TSlrepresented by the 
planar section in which S„ ~ S l,o- Consequently, this allows one to vary the relative contribution to Sc, and thus to 
^.15, of the cusp in which S„ > 5l,o- 

Fig. [T5] shows /0.15, calculated using eqs. (I25]l-(i26]l. as a function of IdSifl for four values of a. At smaller Ic 
close to the width of the reaction zone for a given a (see Fig. [TJb). Iq.\s essentially reflects only the values of S„ in 
the immediate vicinity of the tip of the cusp, where S„ is close to 5*. On the other hand, at larger Ic, the surface 
area of the planar flame region, in which S„ ~ S ifi, is larger, and so this region represents a greater fraction of the 
total surface area of the flame configuration shown in Fig. [12] As a result, the contribution of the planar flame starts 
to dominate that of the cusp, which causes /0.15 — » 1 as seen in Fig. [15] Note that at lower a, the region in which 
S n > S L,o extends progressively further from the cusp tip. In particular, at a = 1°, /0.15 is substantially greater than 
unity even at Ic » 6l.o- At the same time, already at a = 4°, /q 15 drops to < 1.04 at Ic ~ SSifl which shows that S„ 
recovers its laminar value within a few laminar flame widths from the cusp tip. 

These results show that the formation of a cusp due to the collision of planar flame sheets produces values of /0.15 
substantially larger than unity and comparable to those observed in the simulations presented here. As a result of a 
higher local flame speed in a cusp, its contribution to the global turbulent burning speed is disproportionately large 
compared to the fraction of the instantaneous flame surface area in it. A more complex flame configuration, consisting 
of multiple planar flame sheets that approach each other and collide forming multiple cusps, can then be considered. 
In this system, for instance, if the flame surface density is increased, flame sheets will get closer to each other causing 
flame collisions to become more frequent and the resulting cusps to become more numerous. Consequently, the 
fraction of the flame surface area contained in cusps will increase and /0.15 will grow, which is demonstrated by the 
increase in /0.15 with decreasing Ic. 

This simplified model demonstrates the mechanism through which S t can produce an exaggerated response to the 
increase in At- Yet it does not take into account the full complexity of an actual turbulent flame. Such flame does 
not consist of perfectly planar flame sheets that merge, but instead it is constantly wrinkled and folded on a variety 
of scales (Fig. fTTT i. As a result, cusps will form in a multitude of configurations. For instance, in addition to the 
purely 2D situation considered here, there will also exist 3D flame collisions which will result in even higher local 
flame speeds in the cusps and, thus, will have an even larger effect on the magnitude of /0.15. The probability of the 
formation of such more complex flame collisions, however, will rapidly decrease with the increase in their complexity 
is. 23|. Nevertheless, in order to predict a particular value of /0.15 which can be expected in the turbulent flow of a 
specific intensity, it is necessary to understand the types of all cusps which can form, the local flame speed of each 
type, its probability of formation, and, thus, the contribution of each type to the exaggerated response of St- Such 
detailed analysis is the subject for future studies. 



6.4. Effect of cusps on the diffusion velocities in the turbulent flame 

Changes in the local flame structure (Figs. [T2]and[T3ll in cusps will result in modified local diffusion velocities, 
D|VF|. This is the direct consequence of focusing of the diff'usive fluxes which alters the molecular and thermal 
transport in regions of high flame curvature. It is, therefore, instructive to compare the distribution of the diffusion 
velocities in the idealized cusps considered above and in the actual turbulent flame. 

We calculated the surface-averaged normalized diffusion velocity, 

VaiY) kD\^ndA 

Sl AySi ' 

for a set of discrete values of Y in simulation S3 and in the cusp formed at the low inclination angle a = 1° shown in 
the upper panel of Fig. [12] Here the diffusion coefficient D is given by eq. ((7| and the integration is performed over 
an isosurface of a given Y. In S3, the values of Vd{Y) change with time due to the fuel heating by turbulence. In order 
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to take this effect into account, we normalized the instantaneous Vd{Y) by the corresponding actual instantaneous 
laminar flame speed S l, as was discussed in § 12.21 The resulting distributions were subsequently time-averaged. In 
the idealized cusp calculation, S l - S l,q at all times, and the instantaneous flame structure shown in Fi g . [TSlrepresents 
a steady state. The obtained Vd{Y)/S l for both calculations is shown in Fig. [16] We also show two distributions of 
Vd(Y) representing the exact planar laminar flame solutions constructed for two fuel temperatures, namely To - 293 
K and 320 K, which correspond to the average fuel temperatures at the beginning and at the end of S3. These were 
also normalized by the corresponding S l. 

Fig.[T6lshows that the flame configuration in Fig. [12] (upper panel), which consists both of a cusp and the extended 
planar flame regions, has diffusion velocities in the reaction zone < 10% lower than in the planar laminar flame. At 
the same time, a similar decrease of Vd / S l in a laminar flame is produced as a result of fuel heating to 320 K. On the 
other hand, the time-averaged distribution of Vd/S l inside the turbulent flame is also similarly lower in comparison 
with the base laminar case at fuel temperature To - 293 K. This fact, along with the Icr standard deviation of the 
instantaneous values of Vd/S l shown as the shaded gray region, indicates that the observed decrease of the diffusion 
velocities in S3 is consistent with the combined effect of fuel heating and the formation of cusps. 

6.5. Criterion for onset of the regime of flame evolution influenced by cusps 

Properties of cusps discussed above allow us to determine when the evolution of a turbulent flame can be expected 
to become affected by cusp formation, thereby, leading to values of /q 15 substantially larger than unity. The results 
presented in this paper show that such high /0.15 require a large number of cusps to exist or, in other words, the flame 
must have large curvature ~ \I6l over a significant fraction of its total surface. This means that the flame must be 
tightly folded with a characteristic separation 1 /Ej- comparable to the full flame width If, which is indeed the situation 
observed in the simulations discussed here. 

Consider now the two channels of cusp formation discussed in § 16. II and illustrated in Fig.[TT] Nonlocal flame 
collisions cannot create a tightly packed turbulent flame since, by definition, they fold the flame only on large scales 
causing flame sheets to move toward each other. Large curvature forms only at the moments of collision after which 
the large speed of cusp propagation, Uc, causes rapid annealing of the flame surface. Consequently, a tightly packed 
flame can form only if local flame collisions are an efficient process capable of folding the flame on scales ~ 6l not 
only in isolated regions but consistently throughout the flame surface. 

It was discussed in § 16.21 (Fig. [I4l i that Uc, unlike S*, is much greater than Si even at large flame inclination 
angles, and, thus, it increases gradually with curvature. Therefore, Uc is the primary factor responsible for smoothing 
the flame surface while it is being wrinkled by turbulence. In particular, as the turbulent motions U^ shown in Fig.fTTh 
fold the flame thereby progressively increasing its curvature, they must overcome the straightening action of the 
growing Uc- The formation of a cusp, as well as its structure, are then determined by the balance between U,\ and Uc- 

Next consider a section of an idealized curved flame front, schematically shown in Fig. [17] containing a cusp 
formed through a local flame collision (Fig.fTTh'). The flame is perturbed with a wavelength Ac and an amplitude Ic- 
This type of structure was considered by Zel'dovich 1,41.1 as a model of a flame formed under the action of the Landau- 
Darrieus instability in order to analyze the stabihzing effect of cusps. The same structure, however, can emerge under 
the action of any destabilizing process, e.g., the Rayleigh-Taylor instability or, in our case, turbulence. 

In this case, the rate of decrease of the cusp amplitude is |4ll] 

m =-.4-^.1). (28) 

\ dt I- \ sm a I 

The angle a and, thus, (dlcldt)-, will change with the cusp amplitude. In order to relate a and Ic, a parabohc shape of 



the flame was assumed in 14 ill , and eq. (l28T l can then be rewritten as 

(dlc\ „„ /; 



/ die \ C 



The net rate of change of Ic is determined by the balance of the stabilizing process described by eq. ( i29l ) and a 
destabilizing process, namely 

die Ic 

-^^^^-^SL^. (30) 

dt A^ 
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In the turbulent flow, flame perturbations can be assumed to grow linearly in time as they are being stretched by the 
turbulent speed Ua^ characteristic of the scale A^. This gives the growth rate *P = U^^. By setting dlddt - 0, the 
limiting value of IdAc) can be determined. It was shown by Khokhlov I(23i1 in the context of Rayleigh-Taylor-unstable 
flames that nonlinear flame stabilization due to cusp propagation can be expected to fail once Ic > Ac- Therefore, by 
setting Ic = /le, = t/^^., and dlddt = in eq. (l30t . the critical value of the turbulent velocity on a given scale Ac can 



be found 123 1 



f/i, =85z,. (31) 

This is the value of Ua, that is needed to overcome the stabilizing efifect of cusps and, thereby, to deform the flame on 
scale Ac- This result was used to demonstrate that, in order for the flame to be unstable at a wavelength Ac under the 



action of turbulent motions, the turbulent speed on that scale must be substantially larger than S l 112311 . 

This shows that given a specific turbulent cascade, the flame will be stabilized on scales smaller than some critical 
wavelength A"" on which the condition in eq. dSTT l is satisfied. In order for the cusps to have a pronounced effect on 
the flame evolution, the flame must be folded on scales comparable to the full flame width If. Therefore, by setting 
/i"" - If and using eq. OTT l. we find that the flame will be curved on scales ~lf when turbulent velocity on that scale 
becomes ~ 85 l. In a turbulent flow with the Kolmogorov energy spectrum, the corresponding critical value of the 
integral turbulent velocity then is 

* 85l(-J . (32) 

This gives the critical values of the Karlovitz number, /Ta' ''", and the Damkohler number. Da"", written using their 
traditional definitions |i5|] 

^«-'=l£.fiL]''^.8'/^«20, 

^" ,c w,.^/ (33) 



Tf IfUl 

Here Tf = If/Si is the characteristic flame time, t,j is the Kolmogorov time, tt = l/Ui is the characteristic turbulent 
time on the integral scale, and Lc = 1{S l/Ui)^ is the Gibson scale. 

Fig.[T8]presents a traditional combustion regime diagram |5]. The orange region shows the range of the regimes, 
in which /0.15 is expected to be larger than unity, according to eq. ( l33T l. The soUd red square corresponds to the 
simulations discussed here. We also determined the time-averaged value of /0.15 in a simulation similar to S2, but with 
approximately half the turbulent intensity which, however, was still above the Ka'^'" x 20 line. In this calculation, 
represented with the open red square in Fig.[T8] we found /0.15 ~ 1.1, compared with /q 15 » 1.14 found in simulation 
S3 (Table|4|. Details of this calculation will be presented in a separate paper. 

Despite the simplicity of the cusp model used here, eqs. (|32]|-(|33]| provide a rather accurate criterion for the onset 
of the regime of the turbulent flame evolution in which the contribution of cusps is expected to become important. In 
particular, the decrease in the value of /q 15 between the two turbulent intensities considered in Fig. [18] suggests that 
below the Ka'^'" ^ 20 line /0.15, on average, deviates from unity at most by a few percent. 



6.6. Effect of different turbulent intensities and system sizes 

For the given reaction-diff'usion model and fuel properties, the only two free parameters in the system considered 
here are the turbulent intensity and system size represented by the turbulent integral velocity, Ui, and scale, I. Con- 
sequently, the question arises, how will the distributions of St/Sl and, more importantly, of /q 15 as a function of 
Ao.is/L^ shown in Figs.|6j5 and Fig.[8]vary with the change in / and f//? 

The results presented here and in [ 1 1 ] demonstrated the following key property of the flame interaction with high- 
speed turbulence. Turbulent motions are much more efficient at folding the flame with increasingly greater curvature 
than at disrupting its internal structure. This creates tightly packed turbulent flames which, in regions of low curvature, 
retain locally the laminar structure and velocity. Therefore, an increase in Ui creates both larger flame surface area. 
At, and density, Ej-. 

On the other hand, / does not affect Ej-. The analysis in § 16.51 showed that flame folding on a given scale is 
controlled only by the turbulent velocity on that scale rather than by the total range of scales present in the system. 
Instead, a larger system size results in a wider turbulent flame and, thus, in larger values of At. This follows from the 
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fact that, in the high-speed regime, the turbulent flame width, 6t, is primarily determined by / (and, therefore, by the 
driving scale, L) rather than by Ui, as was discussed in [TT] based on the comparison of the results of that work with 
the results of Aspden et al. A similar conclusion was also reached by Peters [5] (cf. eq. (2.175) therein). 

At the same time, a given Ui and / is not characterized by a single St, Aj, and Z^. Instead, a variety of different 
flame configurations are realized in the course of flame evolution, as was seen above (cf. Fig.|6]3). These configurations 
have different Aj, which is the primary factor determining the distribution of St- More importantly, however, at a 
given At, there are also different numbers and types of cusps present in the flame. This creates a scatter in the 
distribution of St/S l vs. At/L^, which in the absence of cusps would lie on the St/S l = AtIL? line. These cusps 
are the result of the combined action of the two cusp-formation channels discussed in § 16.11 Consequently, this scatter 
reflects the dependence of the two types of flame collisions on At and Sj-. 

The efficiency of local flame collisions (Fig. [TTh ) is independent of either At or Sj-, since this process of cusp 
formation is controlled only by the intensity of turbulent motions which stretch and fold the flame at a certain point 
of its surface. As a result, the probability of cusp formation per unit flame surface area in this process does not 
increase with At or 2^. Consequently, if only local flame collisions were present in the system, than the fraction of 
At contained in cusps would grow proportionally to the increase in At, and the resulting /q.is would not change. 

On the other hand, the efficiency of nonlocal flame collisions (Fig. [TTb ) does depend not only on Ui, but also 
on At and S7. To show this, the following qualitative model can be used (see also 123']). In this process, cusps 
are created due to collisions of flame sheets with very low curvature. Therefore, consider an idealized structure 
consisting of planar flame sheets discussed in § 16.31 If the flame surface density of this configuration is S7, then 
the average flame separation is 1 /Sj-. Flame sheets will move toward each other with the speed S l + Ut, where Ut 
is some characteristic turbulent velocity responsible for the advective transport of the flame. Since Ui is the largest 
velocity of coherent turbulent motions, then Ut can be approximated with Ui. Consequently, the flame sheets will 
merge and annihilate in the time tc ~ 1/(2^(5 z, + Ui)). The total number of collisions and, thus, the total number of 
resulting cusps, Nc, will depend on the total surface area of the flame sheets times the frequency of their collisions, i.e., 
Nc oc At/Ic <^ At'^t(S L + Ui). Since At and I.t are wefl correlated (§ l3.3l l. At can be viewed as simply proportional 
to Lt, which give^ A^c (x (S 1 + Ui)A\ oc (S 1 + Ui)1^. Therefore, the number of cusps formed per unit flame surface 
area will be 

^ <x (5l + u)At ^[Sl + Ui)^T- (34) 

These considerations show that at a given / and Ui, the distribution of values of /0.15 has two contributions: one due 
to local flame collisions independent of At and I,t, and the other due to nonlocal flame collisions, which increases 
with At and Ej-. In smaller systems and at lower turbulent intensities, or, equivalently, at lower At and E^, the 
first process is dominant. This can be seen in the distribution of /0.15 vs. Aq.is/L^ in Fig. [8] which shows almost 
no dependence on At- Note, however, a weak increasing trend suggested by the least squares fit. In particular, the 
lowest values of /q 15 become progressively larger with increasing Aq 15 approaching 1.1 as Aq.is/L^ — > 5. While 
the statistical significance of this trend must be confirmed through further studies, this trend would be indicative of 
the contribution of nonlocal collisions. The analysis given above suggests that such dependence of /0.15 on At will 
become pronounced in larger systems or at larger turbulent intensities. 

The fact that /q 15 must depend on At also follows from a different argument. Consider first a limiting value of 
Ej-. The fact that the flame is not an infinitely thin surface means that such a limit does exist. Once E^ becomes large 
enough so that all of the flame surface comes into contact with itself, E?- cannot increase further This maximum value 
E'™' can be estimated by assuming that the minimum flame separation is equal to the full flame width If. Ther0 

ymax _ ymax ^ 

If 



'This expression is very similar to the destruction term often used in the balance equation of the flame surface density (see |44] for a review of 
a number of such models). This is also analogous to the rate of collisions of gas molecules, with Y,t playing the role of number density and the 
turbulent intensity playing the role of temperature in determining the speed with which constituents approach each other. 

'"Note that Z'™' is different from typically used in combustion research and discussed in § 13.21 
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Using the definition of given by eq. ( fT4b , the corresponding hmiting value of the flame surface area then is 

jj^max j^max j j 

T^^^ STcJr = Ci-=(C,C2)-. (36) 

Here, C\ - 6t/L and it typically is w 2 - 4 based on the results of fill The C2 - L/l and in the case of the 
Kolmogorov turbulence C2 ~ 4. In our case Ip ~ 261^ (cf. Fig. 7 in 11 In ). Therefore, based on the values of L and I 
given in Table[3]for the system considered here, S™"' ^1.5 mm"' and A'^^/L^ 7 - 8. 

The situation when all of the flame surface is in contact with itself corresponds to the infinite local flame speed 
and, thus, infinite St- This is equivalent to the complete domination of the nonlocal process of cusp formation when 
the flame sheets collide simply because they are very close to each other and not as a result of the folding action of 
turbulence. Therefore, as At/L^ A"j" jL^, both SjIS l — > °° and /0.15 00. On the other hand, as At/L^ — > 1, 
St/S L — > 1 and, thus, /q 15 1. This shows that /q 15 cannot be constant. Instead, it has to be a monotonically 
increasing function of At, as suggested. 

In summary, the following picture emerges. At a given I, lower turbulent intensities cause the flame to be less 
convolved and, thus, to develop less surface area. As a result, the overall distribution of St/S l vs. At/L- would 
shift toward smaller values of At/L^ and, thus, St/S l. Furthermore, the efficiency of both local and nonlocal flame 
collisions would decrease. As a result, cusps would become infrequent representing only a small fraction of the 
total flame surface. Consequently, the flame almost everywhere would propagate locally with its laminar speed, thus 
causing the scatter in the distribution of St/S l to decrease and values of St/S l to collapse onto the St/S l - At/L'' 
line in the limit of Ui 0. Note, however, that while at Ui < U"" given in eq. dSTT i local flame collisions practically 
cease to occur, individual cusps can still form through nonlocal collisions. Since these depend on At, /0.15 would 
approach unity with decreasing Ui more slowly in larger systems that have greater A /• than in smaller systems. 

At larger Ui, the behavior would be opposite. Not only would the overall distribution of St/S l vs. At/L} shift to 
larger values of At/ and St/S l, but also the efficiency of both local and nonlocal flame collisions would increase, 
causing a much larger scatter of the distribution. Furthermore, at high turbulent intensities, flame configurations with 
very few cusps and, thus, /(> 15 w 1, would become increasingly less likely to occur Therefore, larger Ui will cause 
not only a higher rate of burning due to the increase in At, but it will also periodically lead to a substantially more 
exaggerated response of St to such increase in At- Finally, dependence of /q 15 on At would be more pronounced 
due to the increased contribution of nonlocal collisions and, at a given Ui, larger systems would tend to have higher 
values of /0.15. 



7. Conclusions 

This work continued the analysis of the set of three numerical simulations first presented in fv\^. These calcu- 
lations model the interaction of a premixed flame with high-speed, subsonic, homogeneous, isotropic turbulence in 
an unconfined system, i.e., in the absence of walls and boundaries. The turbulent r.m.s. velocity, t/„„s, is ^ 35 times 
larger than the laminar flame speed, S p. The resulting Damkohler number based on the turbulent integral scale and 
velocity is Da = 0.05. 

It was demonstrated in yjj] that this system represents turbulent combustion in the thin reaction zone regime. Even 
in the presence of such intense turbulence, the flame brush consists of the highly convolved flame with its reaction- 
zone structure virtually identical to that of the planar laminar flame and with the preheat zone broadened by a factor 
X 2. The fact that turbulence is unable to penetrate and disrupt the internal structure of the reaction zone suggested 
that the flame must be propagating locally with the speed S l 1 1 !]• This raised the following question: is the turbulent 
flame speed, St, in the thin reaction zone regime completely determined by the increase in the flame surface area. At, 
as it is the case in the wrinkled and corrugated flamelet regimes? Here we summarize the main findings of our study 
of this issue. 

Analysis of the area and density. Ay and Sy, of the fuel mass-fraction isosurfaces showed that the flame, folded 
inside the flame brush in the presence of high-speed turbulence, cannot be viewed as a thin uniform structure in which 
all isosurfaces are parallel to each other. Diff'erent regions of the flame have quite different response to the action 
of turbulence. In the higher-resolution calculations S2 and S3, both Ay and Ey, on average, increase monotonically 
through the flame with increasing Y, which leads to the distinctive inverted-S shape of the time-averaged distributions 
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of Ay and l^y Consequently, isosurfaces of higher Y are folded by turbulence on progressively smaller scales. This 
causes substantially finer wrinkling of the flame surface on the fuel side than on the product side, as observed in fill]. 
Furthermore, in the reaction zone. Ay grows due to both much tighter folding of isosurfaces and the increase in the 
overall width of the flame brush, with the first process being more dominant. On the other hand, in the preheat zone, 
isosurfaces of Y are packed primarily on small scales without contributing to the increase of dr- 

These properties of the distributions of Ay and Ey showed that, in the thin reaction zone regime, the definition of 
the flame surface area must be revisited before the relation between At and S j can be considered. In particular, it must 
be determined which value of Y fully and accurately characterizes the evolution and global properties of the turbulent 
flame. To answer this question, we analyzed the correlation of Ay with quantities that characterize both the energetics 
and the global structure of the turbulent flame, namely its speed, Sj, and width, Sj- This analysis demonstrated that 
global properties of the turbulent flame in the thin reaction zone regime are accurately represented by the structure 
of the region of peak reaction rate. For the reaction-diffusion model used in this work, this corresponds to Y ^ 0.15. 
Therefore, the isosurface of this value of Y must be viewed as the flame surface and Aj - A0.15. Correspondingly, the 
flame surface density is E7 - S0.15. 

Large values of ILj developed in the course of the system evolution, on one hand, and the absence of any broad- 
ening of the reaction zone observed in 111 ill , on the other, demonstrated the following important fact. High-intensity 
turbulence is much more efficient at tightly packing the flame inside the flame brush than at disrupting and broadening 
its internal structure (also see 15|]). 

The distribution of 5r/5L as a function of A^/L^ and, more specifically, the ratio /q 15 = {Sj IS l)I{At I L?) showed 
that, in the thin reaction zone regime, the magnitude of 5 is controlled primarily by the increase of the flame surface 
area. In particular, on average Iq is ~ 1.14. As a result, a variety of different flame configurations realized over 
time, with both A7 and varying by over a factor of ^ 3, cause St/S l also to vary over a broad range of values 
* 1.6-5.5. Such significant changes in the burning rate demonstrate that, in order to properly characterize the turbulent 
flame evolution in this regime, it is not sufficient to consider only the time-averaged St- Instead, the full PDF of the 
instantaneous values of St must be determined. 

At the same time, the distribution of St/S l exhibited another important characteristic. In the course of the system 
evolution, St/S l varies from within only a few percent of At/ to as high as 30% larger than At/L^- This led to the 
following key conclusion of this work. In the thin reaction zone regime given a sufficiently high turbulent intensity, 
St periodically exhibits a substantially exaggerated response to the increase in Aj-. Such accelerated burning in the 
turbulent flame means that an additional mechanism must augment the effect of the increase of the flame surface area. 

Our analysis showed that tightly packed flame configurations, produced by high-speed turbulence, have the aver- 
age flame separation of only a few full flame widths. This results in frequent flame collisions that lead to the formation 
of regions of high flame curvature > 1/6l, or "cusps." The resulting significant focusing of the thermal flux over an 
extended region of the flame surface can significantly increase the local burning speed, S„,m the cusp over its laminar 
value, S L. These large values of S„ cause the contribution of cusps to the total 5 7 to be disproportionately large com- 
pared to the flame surface area they contain. This provides a natural mechanism to accelerate burning in the turbulent 
flame above what can be attributed to the increase in At, even in reactive mixtures characterized by Le - 1. The 
increase of S„ in cusps is inherently local, and it does not require the flame to be broadened and its speed increased 
by small-scale turbulent transport, in agreement with the results of [11]. 

Therefore, our results indicate that cusp formation is a crucial process in the turbulent flame evolution responsible 
for controlling St- It is capable of accelerating burning in addition to the two processes originally suggested by 
Damkohler [1], namely the increase of At and the enhancement of the local flame speed by turbulence. Furthermore, 
both the criteria given by eqs. (I3ni-(l33ll and the results of numerical simulations show that effects of cusps become 
pronounced well within the thin reaction zone regime and, in particular, at /Ta^"' > 20 (Fig.fTSTl. Above this critical 
value of Ka, increase of St over S l, on average, exceeds by > 10% the corresponding increase of At- This shows 
that a substantially accelerated turbulent burning can arise due to cusps at much lower turbulent intensities than would 
be necessary to disrupt the flame and increase its local speed by small-scale turbulence. 

Above Ka'^'" (Fig.fTSTl. the flame evolution likely remains substantially affected by cusp formation at all turbulent 
intensities, as long as there is folding of the flame by large-scale motions. This is the result of the fact that turbulence 
is more efficient at packing the flame than at broadening it. In particular, at high enough values of Ui, turbulence 
will eventually break the internal flame structure, which will increase S„- This will, in turn, increase the right-hand 
side of eq. (ISTT l. thereby, enhancing the stabilizing effect of cusps. On the other hand, the flame width will also grow. 
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increasing the critical flame separation necessary for the onset of the regime influenced by cusp formation. Since S„is 
determined by small-scale turbulence while flame folding is governed by the faster large-scale motions, it is unlikely 
that Sn can grow fast enough to compensate for the increase in both If and the turbulent speeds which fold the flame. 
Consequently, at higher Ui, i.e., even in the broken reaction zone regime, cusps most likely continue to provide a 
significant enhancement of St- This issue, however, requires further investigation in future studies. 

The effect of flame collisions on the turbulent flame speed suggests that, even at high turbulent intensities, S t does 
not become independent of the laminar flame speed. Analysis given in § l6.5l and. in particular, eqs. (ISTT i- ( l32l l. showed 
that flame folding by turbulence, and thus the efficiency of cusp formation, is controlled by S l. For instance, at the 
same Ui, the reactive mixture with higher S l will develop a less convolved flame with fewer cusps than the one with 
lower S L. Consequently, the exaggerated response of 5 r to the increase of Aj- due to cusps will vary with S l. 

The results presented in this work show that for a large range of turbulent intensities and system sizes, knowledge 
of A 7- is not sufficient to predict the magnitude of St, even if the reaction-zone structure remains unaffected by 
turbulence. For instance, in the regime considered here, using At as a guide would result in errors as high as 30%. 
Such errors are quite substantial, given that the flow evolution is typically very sensitive to the rate of energy release. 
Therefore, it is particularly important to account for the effects of cusps in subgrid-scale models that primarily focus 
on determining the evolution of At (e.g., |23]). Furthermore, when constructing such modified subgrid-scale models, 
flame propagation can no longer be viewed as a local process. In particular, S „ at each point of the flame surface is no 
longer determined only by the local thermodynamic state of the flow or by turbulent motions on scales A < Sl- Instead, 
in the presence of cusps, S„ is also controlled by long-range velocity correlations, which produce flame collisions and 
can span the full size of the system. 

The formation of cusps and the resulting rapid flame propagation in certain regions are a crucial part of the 
turbulent flame-brush evolution in the high-speed regime. Thus, properly capturing this regime in numerical models 
requires the domain size to be larger than the integral scale in order to accommodate the folding of the flame by 
turbulence. Making the domain smaller than / would significantly hamper this process, while making the domain 
smaller than ~(5/, would completely eliminate it. 

The increase of S„ in cusps is discussed here for Le = 1. It is well known, however, that when Le + 1, 5„ varies 
with the flame curvature, even when it is <K 1/^/, |3|. Therefore, due to the large curvature in cusps, any imbalance 
between thermal and diffusion fluxes can significantly exacerbate or suppress the enhancement of S„, depending on 
the value of Le. Therefore, the role and properties of cusps in non-equidiffusive reactive mixtures must be investigated 
in future studies. Furthermore, it is important to extend the results presented in this work by considering the effects 
of detailed chemistry, in particular, to address the role of flame quenching and re-ignition at such high turbulent 
intensities that are necessary for efficient cusp formation. 

Finally, two possibilities exist for the flame evolution in the thin reaction zone regime. If the turbulent intensity is 
not too high, an equilibrium is established between flame-surface creation by turbulence and its rapid destruction in 
cusps. This results in the turbulent flame propagating in a steady state, which is the situation observed here. On the 
other hand, if the turbulent intensity continues to rise, eventually 2^ and At can approach their limiting values given 
by eqs. (l35l l -(l36ll. In this case, both St/S l and St/At can become arbitrarily large, as was discussed in § 16.61 Such 
singular behavior suggests that in reality, unless turbulence substantially alters the flame properties, the steady state 
must cease to exist and the system must undergo a qualitative transformation in order to accommodate the rise in St 
or, equivalently, in the rate of energy release per unit volume. Such qualitative change may indicate the transition from 
the deflagration to a detonation. Detailed discussion of this nonsteady regime will be presented in a separate paper. 
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Table 1. Nomenclature 



Ay = A(Y) 


surface area of the Y isosurface 


Al 


surface area of laminar flame 


At = Ao.15 


surface area of turbulent flame 


B 


pre-exponential factor 


c„ 


specific heat at constant pressure 


D 


molecular ainusion coemcient (eq.[7J 




molecular diffusion constant 




tiirhiilpnt HifFiTiinn popffipipnt Tpn FTli 


Da 


namkhhler number (1 II I's , a1 


I 


ctrf»tr*ri Tur'tcw (f^n 11 n 
bLlCLUil IdL-LtJI V^M-^^ 




ratio of .S T- / "f , to 4 v / / ^ (eci FTfill 


Ka 


Karlovitz number, {If/0 Wi/^l) 


Ka"" 


critical Karlovitz number for cusp formation (eq.[33l 


I 


integral scale 




characteristic cusp length 


If 


full width of laminar flame, ^ 25^ 


r 

1^ 


domain width, energy-injection scale 




Cjibson scale, /(oL.o/t//) 


Le 


Lewis number, 'K/D 


M 


molecular weight 


Map 


Mach number m fuel, (7 (yro/Po) 


Map 


Mach number in product, {/(y/'o/Pf.o) ''^ 


n 


Temperature exponent 


Pq 


initial fuel pressure 


q 


chemical energy release 


Q 


activation energy 


Sl 


actual instantaneous laminar flame speed 


S L.O 


initial laminar flame speed 


Sn 


local flame burning speed 


s: 


maximum local burning speed in the cusp 


St 


turbulent flame speed 


t 


time 




time of ignition 




total simulation time from tj^„ 


To 


initial fuel temperature 




initial post-flame temperature 


U = (Mj-, Uy, U-) 


flow velocity 


u 


turbulent velocity at scale L 


U, 


speed of cusp propagation (eq.l23t 


u, 


integral velocity 


Uf 


critical value of Ui for cusp formation (ea.l32t 


^ nus 


turbulent r.m.s. velocity 




turbulent velocity at scale St o 


Vd 


diffusion velocity 


X = (x, y, z) 


spatial coordinate in the domain 


Y 


fuel mass fraction 


Y 


reaction rate (eq.[8} 




left z-boundary of the flame brush (eq.llOb 




right z-boundary of the flame brush (eg.llOt 
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Table 1 (cont'd) 



Zt,o 




initial flame position along z-axis 


a 




inclination angle between planar flame sheets 


7 




adiabatic index 


Sl 




actual instantaneous thermal width of laminar flame 






initial thermal width of laminar flame 


8j 




turbulent flame width 


Ax 




cell size 


Ktc 




time lag of maximum cross-correlation 


V 




Kolmogorov scale 


s 




energy-injection rate 






thermal conduction constant 


K 




thermal conduction coefficient (eq.|7]l 


Po 




initial fuel density 


Pp.a 




initial post-flame density 


I.y = 




surface density of the Y isosurface (eq.ll4ll 


I.T = 


2o.l5 


surface density of turbulent flame 


Ted 




eddy turnover time, L/U 






time averaging 



Table 2. Input model parameters and resulting computed laminar flame properties 





Input 


To 


293 K 


Po 


1.01 X 10* erg/cm^ 


Po 


8.73 X lO-'' g/cm^ 


7 


1.17 


M 


21 g/mol 


B 


6.85 X 10'2 cmV(g s) 


Q 


46.37 RTo 


1 


43.28 RTo / M 




2.9 X 10-^ g/(s cm K") 


Do 


2.9 X 10-^ g/(s cm K") 


n 


0.7 




Output 


Tp,o 


2135 K 


Pp.0 


1.2 X 10-* g/cm^ 


Sl.q 


0.032 cm 


S L,0 


302 cm/s 
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Table 3. Parameters of simulations^ 



SI S2 S3 

Mesh 64 x 64x 1024 128x 128 x 2048 256 x 256 x 4096 

L 0.259 cm = SSt.o 

Ax 4.05 X 10"^ cm 2.02 x lO^^ cm 1.01 x 10"' cm 

Sl,o/Ax 8 16 32 

Zt,o 1.95 cm = 7.52L 

s 1.26 X lO' erg/(cm3 s) 

Us 4.53 X 10' cm/s = 155 l,o 

U 9.07 X 10^ cm/s = 305 l^o 

Urm, 1.04 X 10* cm/s = 34.485 L,o 

Ui 5.60 X 10^ cm/s = 18.545 l^o 

/ 6.04 X 10-2 cm = 1.87<5l,o 

Ted 2.86 X 10"' s 

tign 3.0Ted S.Ot^j 2.0Ted 

ttotal 16.0Tj,d 

Da 0.05 

Lg 9.47 X 10-« cm = 2.96 x 10-'^Sl,o 

MciF 0.25 

Map 0.09 



Parameters common to all simulations are shown only once in S2 column. 



Table 4. Time-averaged properties of the turbulent flame brush" 





StISl 


0{6tI6l) 




0(5j-/5i,) 


Ao.l5/i' 


0(Ao.i5/L2) So.u, mm-i 


C»(Eo.i5) 


^0.15 


O(/0.15) 


SI 


20.03 




5.19 




3.91 


0.74 




1.31 




82 


17.39 


2.21 


3.91 


2.19 


3.23 


2.84 0.67 


4.07 


1.20 


1.50 


S3 


16.66 




3.55 




3.12 


0.67 




1.14 





^ Time-averaging for all variables is performed over the time interval [iTgd - 16x^^1. 
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Figure 1 Instantaneous kinetic energy spectral density for simulations SI (black), S2 (red), and S3 (green) at a time 
immediately prior to ignition. Shaded regions illustrate scales associated with thermal width 6lx) and full width 26lx) 
of the laminar flame. Vertical dashed lines show wavenumbers corresponding to the Taylor microscales in the product, 
Ap, and fuel, Af, as well as the Kolmogorov scale in the product, rjp, based on the value of the viscosity coefficient 
corresponding to the Sc - Pr - 1 condition. The wavenumber corresponding to the Kolmogorov scale in the fuel, 
T/F = 1.18 X 10 cm, is located outside the range of the graph. (Reproduced from 111 ill .) 
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Figure 2 Isosurfaces of Y in simulation S2 at f = lir^d (cf. Fig. 3, middle row, left panel in iTTj). Isosurface values are 
0.05 (red), 0.6 (green), 0.95 (blue). Red and green isosurfaces bound the reaction zone of the flame. Green and blue 
isosurfaces bound the preheat zone. The zo,,,,,,, and z\,max mark the flame-brush bounds. The zo ,,,,,., and zi^min indicate. 



respectively, the maximum extents of product and fuel penetration into the flame brush (see Appendix A for further 
discussion of zo.max and z\,mm)- (Reproduced from I Till .') 
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y Y 

Figure 4 Evolution of the normalized flame surface area (a) Aq qi/L? and (b) Ao.gg/L^. Evolution of the flame surface 
density (c) Sq.oi and (d) Eogg. (e) Time-averaged normalized isosurface area Ay/L?- (f) Time-averaged isosurface 
density Sy. In panels (e) and (f) time averaging is performed over the time interval [iT^d - 16Tgd], circles represent 
calculated values and solid lines are the Akima spline fits. Shaded gray region shows the distribution of the reaction 
rate, Y, in the exact laminar flame solution normalized by its peak value Ymax = 9.5 x 10^ s In all panels: black 
corresponds to simulation SI, red to S2, and green to S3. 
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Figure 5 Correlation between Sy and Ay /L? in simulation S3 for three values of Y. Graphs exclude data from the time 
interval [0 - 2Ted] during which the turbulent flame develops its equilibrium state. Dashed lines show time-averaged 
values of and Ay/L?-. Solid lines show the least squares fit with its slope given in each panel. 
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Figure 6 Correlation between Ay/L^ and St/S l in simulation S3. Ay is considered for three values of Y: (a, b) 
Y - 0.15, (c, d)Y - 0.5, and {e,f)Y - 0.95. In panels (a), (c), and (e) Ay/L^ is shifted in time to the right by the time 
lag Ate shown in each panel. Graphs in panels (b), (d), and (f) are constructed based on the data shown, respectively, 
in panels (a), (c), and (e) excluding the time [0 - 2Ted] during which the turbulent flame develops its equilibrium state. 
In panel (b) each data point is colored according to the corresponding value of E0.15 with the color scale given in units 
mm"' (cf. Fig.|5^). In panels (b), (d), and (f) dashed lines show time-averaged values of Ay/ L? and S t / S l, solid lines 
show the least squares fit with its slope given in each panel, and dash-dot lines correspond to St/S l -Ay/L^. Shaded 
gray regions in panels (a)-(e) illustrate the exaggerated response of 5 7- to the increase of Ay. 
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Figure 7 Correlation between Ay/L^ and 6t/L in simulation S3 for three values of Y. Graphs exclude data from the 
time interval [0 - iT^d] during which the turbulent flame develops its equilibrium state. Dashed lines show time- 
averaged values of Ay/L^ and 6t/L. Solid lines show the least squares fit with its slope given in each panel. Each 
data point is colored according to the corresponding value of Ey with the common color scale in units mm ' shown in 
panel (a) (cf. Fig.|5]l. 



33 




Figure 8 Instantaneous values of /q 15 as a function of Aq. 15 /L^ in simulation S3. Graph is constructed from data shown 
in Fig.|6j3. Each data point is colored according to the corresponding value of So. 15 with the color scale given in units 
mm Dashed lines show time-averaged values. The solid line gives the least squares fit with its slope given in the 
panel. 




Figure 9 Instantaneous values of /0.15 in simulation S3 as a function of the volume-averaged normalized x-component 
of the kinetic energy inside the flame brush and in the fuel immediately ahead of it. The solid line shows the least 
squares fit with its slope given in the panel. See text for further details. 
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Figure 10 Flame collision and the formation of cusps in the turbulent flame. Shown is the flame-brush structure based 
on the isovolume of Y in simulation S3. Bounding isosurfaces represent Y - 0.05 and Y = 0.95 and the flame brush 
is shown from the product side. Upper panel corresponds to the time t = 11.86r(.^, while the middle and lower panels 
show the flame structure, respectively, O.lr^^ = 3 //s and 0.2Tprf - 6 jus later The thin black line, corresponding to 
Y = 0.6, marks the boundary between the preheat and reaction zones. The thin white line, corresponding to F = 0.2, 
shows the location of the peak reaction rate. Regions A, B, and C show the three main stages of the flame collision 
and the formation of a cusp discussed in § 16.11 Note also two elongated regions of flame collision forming near the 
upper face of the domain. 
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Figure 1 1 Schematic of the two main channels of cusp formation, (a) Local flame collisions. Turbulent motions with 
velocity Ua on the scale A fold the flame, gradually increasing its surface curvature at a specific point until it becomes 
~ 1/<5l. (b) Nonlocal flame collisions. The flame is folded by turbulent motions on a larger scale. The curvature of 
the flame surface remains small at all points until the collision of the two flame sheets causes an abrupt formation 
of two cusps. Resulting cusps propagate in both cases with the phase velocity Uc, which depends on the local flame 
curvature. See text for further details. 




a 



Figure 12 Structure of a cusp formed by the collision of two planar flame sheets. Shown is the distribution of Y for the 
flame inclination angle a - 1° (upper panel) and a - 4° (lower panel). Scale of the panel axes is given in units of 5l,o- 
Away from the cusp tip, the flame propagates in the direction normal to its surface with the laminar flame speed, 5l,o, 
causing the tip to move to the right with the speed Uc. Also shown is the variable length, Ic, of the collision region (see 
text). Thin black line marks the boundary between the reaction and preheat zones, while the thin white line indicates 
the region of peak reaction rate. Note the substantial broadening of the flame, and thus the reaction zone, near the tip 
of the cusp ata- 1° compared to a = 4° (cf. Fig.fTOt. 
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Figure 13 Structure of a cusp formed by the collision of two planar flame sheets for four flame inclination angles, 
a. Shown are distributions of (a) Y and T, as well as (b) the reaction rate, Y, along the symmetry line of the cusp 
(cf. Fig. [T2I 1. Dashed lines in both panels indicate the exact planar laminar flame solution. Profiles of T and Y are 
normalized by their respective maximum values in the exact laminar solution, i.e., Tpfl (see Table |2]l and Y„ax - 
9.5xl0'*s-i. 
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Figure 14 (a) Dependence of the normalized speed of cusp propagation, Uc/S l,o, on the flame inclination angle, a. 
Solid line corresponds to the analytic expression given in eq. ( |23] ). red circles show the computed values. Horizontal 
dashed line indicates sound speed in cold fuel, (b) Dependence of the maximum normalized local flame speed in the 
cusp, S*„/S L,o, given by eq. (|24| ). on a. Red circles show the computed values, solid line is the Akima spline fit. 
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Figure 15 Dependence of /0.15, given by eq. ( |26] |. on the normalized length of the flame-collision region, IdSifi, for 
four values of a. The shaded gray area indicates the full laminar flame width If « ISifl- See text for further details. 




Figure 16 Time-averaged distributions of the normalized diffusion velocity in simulation S3 (black line) and in the 
idealized 2D cusp at a = 1° shown in the upper panel of Fig. [12] (green line). Circles are calculated values, solid lines 
are the Akima spline fits. The shaded gray region is the Icr standard deviation of the instantaneous values of Vd/S l 
in S3. Dashed red lines are the exact laminar flame solutions corresponding to fuel temperatures Tq - 293 K (upper 
line) and 320 K (lower line). The shaded orange region shows the distribution of the reaction rate, Y, in the exact 
laminar flame solution (fuel temperature Tq) normalized by its peak value Y,nax - 9.5 x lO'* s See text for further 
details. 
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Figure 17 Illustration of the idealized perturbed flame stabilized by the propagation of a cusp with the speed Uc (cf. 
Figs.[n}iand[12li. 
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Figure 18 Combustion regime diagram according to fs*]. The orange region above the Ka"" « 20 line shows, ac- 
cording to the criterion discussed in § 16.51 the range of the regimes in which the formation of cusps is expected to 
substantially affect the turbulent flame speed. The solid red square corresponds to the simulation S3 presented in 
this work that gives /q 15 a; 1.14. The open red square shows the regime with one half the turbulent intensity in 
which the value of /0.15 ~ 1.1 was determined. Flamelets are typically believed to exist in the regimes below the 
Ka^ - /Ta/lOO = 1 line (Ist]. The traditional form of the diagram was also modified by adding the line Map - 1 
indicating the region of supersonic turbulence in the cold H2-air fuel under the atmospheric conditions. See text for 
further details. 
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Appendix A. Definition of ttie flame surface density 



In the definition of Ey, given in eq. (fT4b . the isosurface area, Ay, was normaUzed by the total volume of the flame 
brush, djL?- Individual isosurfaces, however, do not occupy all of this volume. For instance, this is the case for 
isosurfaces of F = 0.05 and Y = 0.95, which are used to specify St in eqs. (l9ll-(fT0t. simply because, by definition, 
they are separated at least by the full flame width. If, or a larger distance if the flame becomes broadened. In our case, 
Ip » 26l (see Fig. 7 in 1 1 1 ]) is a substantial fraction of 6t- Consider Fig. |2] It can be seen that neither the Y - 0.05, 
nor the Y = 0.95 isosurface extend over the full width of the flame brush, and they indeed occupy a smaller volume. 

In order to consider this more quantitatively, we define zo,max and z\,mm by analogy with zo,mm and zi,max (ea.fTOll. 

= min(z) : Y(x, y, z) > 0.05 V (x, y,z> Zo,max), 
= max(z) : y(x,3;,z) < 0.95 V (jc,y,z < zi,™„). 

In other words, zo^max is the z-coordinate of the rightmost cell with pure product, while z\,mm, respectively, is the z- 
coordinate of the leftmost cell with pure fuel. Thereby, they effectively measure the furthest extent of product and fuel 
penetration into the flame brush. Given the definitions of zo.mm and z\,max (eq.fTOll. zo^min and zo.ma.i bound the volume 
confining the Y - 0.05 isosurface, while zi_„„„ and zi^nmx bound the Y - 0.95 isosurface. If zi^nin < zo,max, these two 
variables bound the region of macroscopic mixing of product and fuel in the flame brush, i.e., the region in which both 
pure fuel and pure product can be found. On the other hand, values zu,max ^ Zi.mm exist only in a weakly wrinkled 
(nearly planar) flame. Fig.|2]provides the illustration of all four of these quantities. 
Next, we define the position of the turbulent flame brush as 

Zt = ^ ■ (A-2) 

Modified zo^max and zi,mm are then defined as offsets of zo,max and zi,mm with respect to zt normaUzed by 6t, namely 

_ ZO,max ~ Zt ^ Zl^min ~ Zt , . 

ZO,max ~ Z ? Zl^min — Z ■ v-^--^/ 

Ot Ot 

Thus, zo.ma.v and zi,„„„ are the relative measure of the extent to which pure fuel and pure product penetrate into the 
flame brush. It follows from the definition (IA.3b . that both quantities take on values in the interval [-0.5,0.5]. For 
example, in a planar laminar flame, they are constant with zo,max = -0.5, 2i,,„!n = 0.5. In the turbulent flame brush, 
Zo,max = 2i,mm = would indicate that both fuel and product reach the midpoint of the flame brush, they are confined to 
the left and right halves of the brush, and there is no macroscopic mixing of fuel and product. This would correspond 
to a weakly wrinkled flame. If zo^max ~ 0.5 and z\,mm ~ -0.5, then both pure product and pure fuel can be found 
throughout the entire volume of the flame brush. 

Fig. lB.19l shows the evolution of zo,max and li_,„;„ for simulations S1-S3. In all cases, both parameters exhibit fairly 
similar behavior oscillating around zero. As the system evolves, it undergoes recurring transitions between periods of 
enhanced fuel-product mixing and episodes of their near complete separation when the turbulent flame becomes more 
planar The correlation between these quantities and 5 r is less prominent compared, for instance, with that between 
Ay and St, although it is possible to associate some peaks and troughs with the corresponding changes in 5 7-. 

It follows from Fig. lB.19l that throughout the course of the simulation, both the Y = 0.05 and Y - 0.95 isosurfaces 
indeed occupy regions smaller than the full flame-brush volume. This is also the case for other isosurfaces (cf. Fig.|2]l. 
Consequently, an argument can be made that when calculating Zy, normalization should be performed not over the 
full volume of the flame brush, but rather over the respective volume bounding a given isosurface. For instance, 
for Y - 0.05, this volume is (zo.max - z%mm)L^ — izo,max + Q.5)6tL^, which would introduce an additional factor 
i/izo,max + 0.5) in eq. (fT4l i. Such definition of Ey would be a more accurate measure of how tightly a given isosurface 
is folded. Fig. lB.19] shows that, on average, zo.max ~ and, thus, the Y = 0.05 isosurface is confined to about half of the 
total volume of the flame brush. Therefore, such modified definition would result in approximately twice higher values 
of Ey. Ultimately, however, we are interested in determining how the overall flame, rather than individual isosurfaces, 
is folded inside the flame brush since only the flame as a whole has actual physical significance. Consequently, the 
uniform normalization over the total volume occupied by the flame appears to be a more physically grounded choice. 

Further motivation for our choice of normalization in the definition (fT4l l of Ey comes when we consider how zo,max 
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and zi,min would change for larger system sizes. Since the Y - 0.05 and Y - 0.95 isosurfaces are always separated at 
least by the full width of the flame. If, then 



Z\,max - Z0.max ^ If, Zl,mm - Zo.min ^ If- (A.4) 

Using eqs. iA.li . and ( IA.3I ). and the fact that in our system If x 26l, inequalities (IA.4I) can be transformed into 
the following conditions which must be satisfied at each moment in time, 

Zo,,„a.x < 0.5 - — , zi,,„m > -0.5 + — . (A.5) 

Ot Ot 

Horizontal dashed lines in Fig. |B.19l show the average limiting zo,max and zi_,„,„ based on the values of 6t/Sl listed in 
Table|4]for each simulation. In particular, zo,ma.x must be <0.4 while 2i,mi„ must be > -0.4, which closely agrees with 
the data shown in Fig. lB.19l 

It was discussed in § 16.61 that 6t increases with the turbulent integral scale, or, equivalently, with the system 
size and energy injection scale ifTll Isll. It then follows from eq. iA.5\ that as 6tISl then limiting values 

of z%max — > 0.5 and z\,min — * -0.5. Therefore, in larger systems, as 6t becomes large in comparison with 6l, the 
volume bounding each isosurface becomes well approximated by the total volume of the flame brush. Consequently, 
in the limit of large values of St, the uniform normalization for all values of Y in eq. (fT4l i becomes equivalent to the 
normalization by the actual volume bounding a given isosurface. 



Appendix B. Induction time for the planar laminar flame 

Fig. lB.20l shows the distribution of induction times t,„j, calculated using the procedure described in § 14.11 through- 
out the planar laminar flame. The solid red line corresponds to the fuel temperature Tq - 293 K and density po (see 
Table |2]l, while the dashed red line was obtained for the higher fuel temperature 320 K and the same density. The 
second temperature reflects the efifect of fuel heating by turbulence in the course of the simulation (§ 12. 2b . 

For comparison, we also show the adiabatic induction time, t"^^^, calculated using eq. ( l20l i based on the exact 
planar laminar flame solutions. Solid and dashed line correspond, respectively, to the same fuel temperatures as for 
Tind- Overall, throughout the reaction zone, i.e., for Y - 0.2 - 0.6, t,„^ and tJ|^^ differ by less than a factor of x 2. 
Moreover, near the peak reaction rate, i.e., at F » 0.2, and at the boundary between the reaction and preheat zones, 
i.e., at y X! 0.55, Ti„d and t"^^ become virtually equal. In the preheat zone, however, they begin to diverge substantially. 
This is to be expected, since t^^^ considers only self-heating of a fluid element by chemical reactions and neglects the 
diffusive transport of heat, which is dominant in the preheat zone. 

Despite the reasonable agreement in the reaction zone between t,„^ and t".^^^, t"^^^ should not be viewed as a model 
of Tjnd, but rather as a useful approximation. In particular. Fig. IB.20l shows that t^^^ does not capture the correct shape 
of the distribution of r,,,,/. The t;,,^/ increases more slowly at larger Y than t"^^^ since slower heating of a fluid element 
by chemical reactions is partially compensated by the larger influx of heat due to the thermal conduction. Since t^J^^ 
does not account for this, its rate of growth with Y is larger in the outer regions of the reaction zone. At the same time, 
since energy balance of a fluid element in the reaction zone is still dominated by self-heating rather than by thermal 
transport, t"^^ remains close to the actual value of the induction time given by Ti„d- 

We also show in Fig. IB.20l the characteristic reaction time, Treac, defined based on the reaction rate (eq. [8) and 
neglecting the change in the fuel concentration. 



'^reaci^ ) — 




(B.l) 



Similarly to t^J^^, T' and p' for a given Y' are determined from the exact planar laminar flame solutions corresponding 
to the fuel temperatures Tq - 293 K (solid line) and 320 K (dashed line) and density po. Unlike t^^^, which is typically 
applicable in the context of an autoignition process, this quantity is often used to approximate the induction time in 
the flame. Nonetheless, for all Y > 0.15, Trmc provides a less accurate approximation of T,„rf than t^^^. 

Finally, blue circles show the correlation time lags Ate calculated in § 14.11 There is close agreement between Ate 
and both t,„,/ and r"^^. Furthermore, when fuel heating is taken into account (dashed lines), all three quantities become 
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Figure B.19 Evolution of the normalized maximum extents of product and fuel penetration into the flame brush zo,ma.t 
and z\,min (see Appendix A for the definition and Fig.|2]for the illustration) for simulations SI (a), S2 (b), and S3 (c). 
The zo,/iifl.v is shown with red line, zi.m,„ - with blue line. Shaded regions mark the extent of macroscopic mixing of 
pure fuel and product inside the flame brush. Horizontal dashed lines show the average limiting values of zo.max and 
Zi.rnm given by the condition (IA.5b and calculated using the values of ^r/^L Usted in Table H] 



even closer. 
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Figure B.20 Distributions of the actual induction times, t,„j, adiabatic induction times, t^^^^, and the characteristic 
reaction times, Treac, in the planar laminar flame. Also are shown the correlation time lags Af^ (cf. Fig. |6l). Red 
circles show the calculated values of Tind, red lines are the Akima spline fits. Green and black lines represent eqs. ( l20l l 
and dB.ll ). Solid lines correspond to the fuel temperature Tq - 293 K, dashed lines to 320 K. Shaded gray region 
shows the distribution of the reaction rate, Y, in the exact laminar flame solution normalized by its peak value Y,„ax - 
9.5 X lO'* s The vertical dashed line indicates the region of peak reaction rate. The horizontal dashed line indicates 
equality between the induction time and the large-scale eddy turnover time, Ted, during which turbulence completely 
reorganizes the turbulent flame structure. 
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